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Mario Quintana - A Cor do Invisi'vel 



Abstract 



This thesis consists of two parts that are independent of each other. 

In the first part I report on a joint work with V. Sidoravicius concerning the 
Activated Random Walk Model. This is a conservative particle system on the 
lattice, with a Markovian continuous-time evolution. Active particles perform 
random walks without interaction, and they may as well change their state to 
passive, then stopping to jump. When particles of both types occupy the same 
site, they all become active. This model exhibits phase transition in the sense 
that for low initial densities the system locally fixates and for high densities it 
keeps active. Though extensively studied in the physics literature, the matter of 
giving a mathematical proof of such phase transition remained as an open problem 
for several years. In this work we identify some variables that are sufficient to 
characterize fixation and at the same time are stochastically monotone in the 
model's parameters. We employ an explicit graphical representation in order 
to obtain monotonicity. This representation has a very useful commutativity 
property that allows direct, constructive approaches. With this method we prove 
that there is a unique phase transition for the one-dimensional finite-range random 
walk. 

In the second part I report on a joint work with V. Sidoravicius, D. Surgailis, and 
M. E. Vares about the Broken Line Process. In this work we introduce the broken 
line process and derive some of its properties. Its discrete version is presented 
first and a natural generalization to the continuum is then proposed and studied. 
The broken lines are related to the Young diagram and the Hammersley process 
and are useful for computing last passage percolation values and finding maximal 
oriented paths. For a class of passage time distributions there is a family of 
boundary conditions that make the process stationary and reversible. For such 
distributions there is a law of large numbers and the process extends to the infinite 
lattice. One application is a simple proof of the explicit law of large numbers for 
last passage percolation with exponential and geometric distributions. 
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Introduction 



In Chapter [T] we study a system of activated random walks ( ARW) on the infinite 
one-dimensional lattice Z. We assume that there are infinitely many particles in 
the system, each of which can be in one of two states: B (active) or A (passive 
or sleeping). Each 5-particle performs an independent, continuous-time, simple 
random walk on Z, with the same jump rate, which we assume, without loss 
of generality, to be equal to 1. We will be treating both the symmetric and 
asymmetric walk cases. When a 5-particle jumps to a site with one or more A- 
particles, any such particle at this site is immediately activated (i.e., switches to 
state B). Each isolated £>-particle goes to sleep (switches to state A), at a rate 
A > 0. From this rule it follows that if two or more particles occupy the same site, 
then they are all of type B or all of type A (the latter situation can only arise 
in the initial condition). According to this rule, at most one S-particle can go to 
sleep per site, and, if undisturbed, remain in state A forever after. 

This model belongs to a broad class of interacting particle systems with conserva- 
tion, which have attracted great interest in physics, probability, and allied fields, 
in part because they afford simple examples of phase transitions in systems main- 
tained far from equilibrium. In this class of models the particles exist in two states 
that may be termed active and passive. Such that activation of a passive parti- 
cle requires the intervention of one or more active ones. This class includes the 
so-called conserved lattice gases [Liib02a| ILub02b| ILH03al ILH03b| IRPSV00] and 



stochastic sandpile models [DTdO02l IDic081 IDAM+OH IMan90[ IMan91] . Such 
models exhibit self-organized criticality [BTW87, IBTW881 IDha99l IGri95] when 
coupled to a suitable control mechanism |Dic02t [DMVZ00J . 

In fact, the stochastic conserved sandpile, generally known as Manna's 
Model |Man90| IMan91j . served as the primary motivation for the present study. 
In infinite volume, this model is defined as follows. Initially there are infinitely 
many particles, distributed in such a way that at each site of Z d we have a Poisson 
number of particles with mean fi > 0. Each site is equipped with an exponential 
rate 1 clock, and each time a clock rings at a site bearing 2d or more particles, 
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2d particles jump to randomly chosen nearest neighbors. Differently from the de- 
terministic Bak-Tang-Wiesenfeld sandpile model [BTW87, BTW88], each particle 
chooses its direction among the 2d possibilities with probability ^, independent 
of any other particle. In contrast to the deterministic sandpile |Dha99] , very lit- 
tle is known rigorously about this system, and the ARW model is a reasonable 
caricature that seems to capture some essential aspects of Manna's model. 

The ARW model may also be viewed as a special case of a diffusive epidemic 
process. In this process, an infected particle performs a simple symmetric ran- 
dom walk with jump rate Db, and recuperates at a given rate, while a healthy 
particle performs a simple symmetric random walk with jump rate Da', healthy 
particles are infected on contact with infected ones. The ARW model corresponds 
to Da = 0. The model was introduced to probabilistic community (in the case 
Da = Db) in the late 1970's by F. Spitzer, but due to its tremendous technical 
difficulties and complexity, remained unsolved until recently, when it was stud- 
ied in detail in [KS031 [KS05j IKS06[ IKS08| . The diffusive epidemic process has 
also been studied via renormalization group techniques and numerical simulation 
[dFLdSHOOl IFMLOlal IFMLOlbl iJanOTl IKSS89L IUvWLH89llvWUH98j . A general 
conclusion from these studies is that there are three distinct regimes of critical 
behavior, for Da < Db, Da = Db and Da > Db- It is not yet clear whether 
the ARW model falls within the remit of the first regime, or, alternatively, that 
Da = marks a special case. 

Here we suppose that each site x G Z initially contains a certain number T]q(x) 
of -B-particles and no A-particles. The t]q{x) are i.i.d., mean //, Poisson random 
variables. 

We are mainly interested in the phenomenon of fixation. By fixation we mean 
that almost surely, for any finite volume A, there is a finite time t\ such that after 
this time there are no 5-particles within A. When there is no fixation we expect 
that there is a limiting density of active particles for long times. 

Numerical analysis and some general theoretical arguments suggest that the ARW 
model exhibits a phase transition in the parameters A and /i, and that there should 
be two distinct regimes: 

Low particle density. There is a phase transition in A in this case. Namely, there 
is < A c < oo such that the system fixates for A > A c and does not fixate for 
A < A c . The threshold A c should vanish when \i tends to zero and it should diverge 
as n approaches the high-density regime. 

High particle density. In this case there is no phase transition and the system 
does not fixate for any value of A. 
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The second regime was studied in }KS06j . where it was shown that indeed if the 
initial density is large enough there is no phase transition. 

On the other hand, the behavior predicted for low particle densities, in spite of 
being nearly obvious, resisted many attempts to be established rigorously. 

In this work we estabilish part of the predictions for in the one-dimensional case. 
The main result is that A c < oo for \i < 1 and A c — > as fi I 0. 

We also prove that A c = oo for /j, ^ 1, characterizing what ranges of \x pertain 
to each regime. Yet two conjectures remain open: that A c > for all fi > and 
A c — > oo as /j, 1 1. 



In Chapter [5] we study a continuous variant of so-called Broken Line Process. The 
meaning in which we use word "continuous" requires more detailed explanation, 
which will be given below. 

Speaking in more general terms, this process might be viewed as a generaliza- 
tion of well known Hammersley Process, introduced by Aldous and Diaconis, and 
independently by Rost. 

The Hammersley Process belongs to the broad class of Polygonal Markov Fields, 
which have been promoted by A. N. Kolmogorov in his later years. However, 
due to technical difficulties arising from the interplay between stochastic geom- 
etry and an intention to construct and preserve a spatial analog of the Markov 
property, this topic became infamous for not only being difficult technically, but 
also conceptually. 

As a consequence, until now it remains poorly understood. A few fundamental 
works were produced by Arak and Surgailis |AS89al IAS89bl IAS89d IAS89dl IA"590l 
IAS91] , which had fair success among the community of statisticians and statistical 
physicists studying random patterns. We find in writings of J. Hammersley some 
indication that he had thought about particular cases of such fields, however 
without developing much of research in this direction. 

Later on, in the context of applications to the First and Last Passage Percola- 
tion, the process which we call nowadays Hammersley Process, was introduced 
by Aldous and Diaconis [AD95] . and Rost [Rosj . and on more general basis stud- 
ied by Sidoravicius, Surgailis, and Vares [SSV99J. Results of |SSV99j were ex- 
tended to the last passage percolation model on the square lattice with geometric 
weights [SEYj. 

In this work we generalize and extend |SSVj to the case of lattice models with 
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exponential weights (and in this sense we use the word "continuous"). We derive 
basic properties of this random field, first of all its existence. Finally we identify all 
possible stationary broken line processes related to both geometric or exponential 
weights to be given by a one-parameter family. As a consequence we obtain the 
law of large numbers for the asymptotic velocity of last passage percolation by 
relatively simple and transparent arguments - results that were previously known 
and are summarized in |JohOO] . Fluctuations, however, are still beyond of the 
reach for these geometric techniques. 



Chapter 1 



Phase Transition for Activated 
Random Walk Models 

The model of activated random walks evolves as follows. On Z d , at time zero 
there is an i.i.d. number of active particles whose expectation is ji. Each particle 
performs a continuous-time random walk with rate one. When a particle is found 
alone at some site, it will change its state from active to passive when an expo- 
nential clock of rate A rings. Once a particle is passive, it no longer jumps. The 
particle becomes active again when some other particle jumps into the same site, 
if that ever happens. 

The phenomenon we consider is that of local fixation. By local fixation we mean 
that at any finite box a.s. there will be only passive particles or empty sites for 
large enough times. Absence of fixation means that at every site there are at least 
two particles for arbitrarily large times. A natural conjecture is that for low \i the 
system locally fixates and for large \i there is no fixation. 

The first difficulty for the study of this model lies on the nature of interactions, 
that even though sound very simple, are quite severe to analyze and forbid the 
use of the most common available tools. Namely the evolution of this model does 
not preserve any type of domination, or, in the language of particle systems, this 
model is not attractive. We identify an observable that suffices to characterize 
fixation and we introduce an explicit graphical representation that recovers the 
monotonicity of this variable. Such representation has also a very useful commu- 
tativity property that strongly favors direct, constructive approaches. With this 
method we prove that there is a unique phase transition for the one-dimensional 
finite-range random walk. 

This chapter is based on a joint work with V. Sidoravicius. 
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CHAPTER 1. PHASE TRANSITION FOR ARW MODELS 



1.1 Introduction 

We study a particle system that is informally described as follows. The system 
starts with an i.i.d. Poisson(yu) number of active particles at each site of the lattice 
S = Z d . As time passes these particles perform continuous-time random walks. 
The particles do not interact except that when there is only one particle at a site, 
this particle may become passive, or fall asleep, at rate A > 0. Once the particle 
is passive it will stop moving until some other particle jumps into that site. We 
say that particles which are active are particles of type B and particles which are 
passive are of type A. In the dynamics the -B-particles jump at rate Db = 1 and 
A-particles jump at rate Da = 0; a .B-particle becomes type A at rate A and an 
A-particle becomes type B immediately if there is some other .B-particle at the 
same site. 

The state of the system is described by an element of S = (Z + x and 
is denoted r\ = (j) (x), V B ( x )) xe g- In this setting r)f(x) denotes the number of 
A-particles (that is, particles which are passive) at position x at time t; r]f (x) 
denotes the number of -B-particles; and rj AB (x) = r]^(x) + rj B (x) denotes the total 
number of particles. 

It is known that there is an underlying probability space where a cadlag strong 
Markov process taking values on an appropriate subset E'CE and corresponding 
to the dynamics described above is defined. This technical issue is discussed in 
Section 11.21 

A natural question that arises is whether the system fixates, which amounts to 
say that almost surely r]f(x) = for all x in a given finite box A and for all t large 
enough. In particular one wonders if there is phase transition, which means that 
there is fixation when the initial density is small enough and there is no fixation 
when it is large enough. More precisely, phase transition means that for each fixed 
A there exists a /x c G (0, oo) such that the systems fixates for \i < \i c and does not 
fixate for \i > fi c . 

A severe difficulty one has to plunge into when approaching this question lies 
on the nature of the interactions that take place. To be more precise, most 
attempts to attack this problem using modern tools and techniques have failed 
for the following reasons. Good estimates that locally solve the problem do not 
suffice; once the evolution halts at a finite region, there is no guarantee that in 
the future more particles from outside will not enter the controlled region. This 
restarts the process and invalidates the argument. A renormalization approach 
is the standard technique to overcome the shortcomings with the local argument. 
Typically one gets better and better estimates for all scales, and finally combines 
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them to obtain a global result from local estimates. This type of approach breaks 
down because they rely on attractiveness, a property that our model unfortunately 
lacks. (Another consequence of non-attractiveness is that it is not even entirely 
obvious that there is only one transition \i c e [0, oo].) 

This chapter is mainly devoted to the proof of 

Theorem 1.1. Consider the activated random walk model on the lattice /L d with 
fixed parameter A > 0. Then there exists ^ \l c ^ oo such that the system fixates 
for ji < \i c and does not fixate for \i > \i c . 

If d = 1 and the random walk consists of finite-range jumps then < fi c ^ 1. 
Moreover, if the jumps are only to nearest neighbors then fi c is 

To overcome the difficulties that arise from the lack of attractiveness, we consider 
an explicit graphical construction that has strong commutativity and monotonic- 
ity properties. This construction will form the core of the proof. 

The fact that fi c < oo is intuitively obvious, since for fi > 1 there are 'more par- 
ticles than sites'. Indeed this proof is quite simple once the framework discussed 
above is set up. Less obvious is the conjecture that \i c should be strictly less than 
one. Another open question is whether there is fixation at fi = /i c . For the totally 
asymmetric, nearest-neighbor walk, a negative answer was given by Hoffman and 
Sidoravicius [HSJ. 

Section 11.21 concerns technical questions such as the existence of the process and 
the approximation of the probability of local events on the original model by a 
large, finite one. We present the graphical construction in Section [T73l and study its 
properties. In Section 11.41 we prove Theorem 11.11 using the background developed 
in the previous sections. Finally, in Section 11.51 we make some general comments 
on generalizations of the proofs presented here as well as some open questions. 

1.2 Existence of the process and approximation 
by finite boxes 

In this section we discuss the construction of a strong Markov process with cadlag 
trajectories whose dynamics correspond to the infinitesimal description given in 
the Introduction. 

On the lattice S = Z d , consider a translation-invariant random walk distribution 
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p(x, y) = p(0, y — x),x,y G S. The formal generator of our process is given by 

(Lf)(v) = E ^i(v B (x))[f(v {x) ) - f(v)] + £ ^ B ( X M X ' y)[f(v {xy) ) -f(v)) (i.i) 



for / : £ — > R, where 



77(a;) , Z X 

(1,0), z = x 



and 



7/(2;), z^x,y 



^ ) («) = < (0,^(a;)-l). 



(0,^(y)+77 B (y) + l), z = y 

Of course the evolution is clearly well defined for any initial configuration in 
E" = {j; 6 E : ^2 X V AB ( X ) < °°}- That is, given a probability distribution P 
of the initial configuration 77 G £", there is a unique Markov process defined for 
£ G [0, 00) with values in E" and having (II. ip as generator. Let P denote the law 
of the process with initial distribution P. 

Given M G N, take P^ f as a product measure on E having marginals P^ I (^r] A (x) = 
0) = 1 V x G S and 



P^(rA*)=n) 



HI < M 



£ (n), \\x\\ > M 



As above, P^ will denote the law of the process with initial distribution Pj^. 

The dynamics we want to construct cannot be defined on the whole space E 
and we ought to restrict the set of allowed configurations. Define a(x) = 
Er=o 2 "> (n) (^0) and |M| = E x& s{\v A (x)\ + \r] B (x)\)a(x). Let £' = {77 G E : 
II77H < 00}. Then E" is dense in E', P M (£') = 1 and straightforward adaptations 
of the Andjel's construction [And82j imply the existence of a unique measure P^ 1 
on Q = D([0, 00), £') with the property that 

W(A) = lim F^(A) (1.2) 

M— >oo 

for any event A such that, for some finite ACS' and < t < 00, A is measurable 
with respect to (v?' B (x)) xeAjSem ■ 
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1.3 A graphical construction that preserves 
monotonicity 

Since the model is not attractive, some type of comparison between different 
parameters (initial density, jump rate, sleep rate), as well as the obtainment of 
bounds, are highly desirable for its study In this section we shall give a graphical 
construction of the dynamics. This construction recovers the monotonicity of 
some variables of interest, allowing the desired comparisons. 

It is simpler, clearer and more convenient to state and prove such properties for 
initial states with only finitely many particles on the lattice. In particular we shall 
construct explicitly. Finally in Section [PI we shall use (jl.2p . together with 
the results presented here, to study the phenomenon of fixation. 

The construction holds on a general graph S and with a general random walk 
distribution. For the sake of clarity we assume S = Z d , d — 1, and consider 
nearest-neighbor random walks with probability p of jumping to the right and 
q = 1 — p of jumping to the left. 

The dynamics is described in a informal way as follows. Every particle waits an 
exponential time of rate 1 + A and when the clock rings this particle will perform 
an action, which may be a jump or an attempt to sleep. But the particle will not 
toss a coin to decide what to do. Instead, it is the site where the particle is at 
that moment that will toss the coin. It seems at first that it makes no difference 
but it is with that intuition that we conceive the explicit construction. In this 
setting we place envelopes at the sites. 

Suppose there are finitely many particles on the system. Let there be a universal 
clock that will ring with the appropriate rate and let there be an i.i.d. sequence of 
labels, independent of the clock. Also, at each site, let there be an i.i.d. sequence 
of envelopes, each one containing some instruction to be performed. When the 
clock rings for the first time, it will ring for the particle indicated by the first 
element in the label sequence and at that moment this particle will perform some 
action. If the particle is passive, nothing happens. If the particle is active, it will 
open the first envelope at that site, burn the envelope and perform the action 
written inside. The instruction may be 'jump to the left', 'jump to the right', 
or 'try to sleep'. So we say there are two types of envelopes: jump envelops and 
sleep envelopes. The envelope is burned regardless if there are other particles on 
the same site and the particle unsuccessfully attempts to sleep. 

This graphical construction, that is in particular a natural coupling, has two nice 
properties that we describe below. 
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The first property is commutativity and says the following. Suppose for a given 
realization of the process (universal clock, label sequence, envelopes) the system 
stabilizes, that is, all particles are passive for large enough time (of course start- 
ing with finitely many particles this happens a.s.). Then by changing the label 
sequence and the universal clock the system will stabilize at exactly the same 
state, except that some particles may be permuted. Furthermore, the amount of 
envelopes burned at each site is also preserved. In Figure [T7TI we show an example 
of a setup given by initial position of the particles and a stack of envelopes at 
each site. The reader may evolve the system as she likes, and the evolution will 
terminate after burning exactly the envelopes indicated. 




Figure 1.1: commutativity. The one-dimensional lattice with the particles' starting 
positions and the envelope sequences. Particles are represented by black balls below 
the axis and envelopes are represented by arrows (jump envelopes) or crosses (sleep 
envelopes) just above the axis. No matter in which order the evolution is performed, it 
will end up burning the envelopes below the bold line. 



So the final state of the system is determined by the initial conditions (positions 
and types of the particles) and the by sequences of envelopes. The the second 
property, monotonicity, states the following. Suppose for some realization of the 
envelopes and initial conditions there is fixation, and take a new configuration 
by deleting some particles on the initial condition, changing the type of some 
particles from active to passive or inserting some sleep envelopes at some sites' 
sequences. Then for this new configuration the system also stabilizes and the final 
number of envelopes that are burned at each site (not counting the ones inserted) 
does not increase. This is illustrated in Figure IT72[ 

Graphical constructions in general have become standard tools in Probability and 
are usually attributed to Harris. Actually it was known before but we have not 
determined exactly when it first appeared. See Harris |Har78] and references 
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x 




Figure 1.2: monotonicity. First, inserting a couple of sleep envelopes has decreased the 
final profile of burned envelopes. Second, starting with some particles passive instead of 
active (squares instead of balls) has caused an analogous effect. Finally, even starting 
with quite a few extra sleep envelopes and some passive particles, the profile does not 
necessarily decrease. 



therein to backtrack this history. In the construction described above only the 
instructions are left at each site; the instant when each step should be performed 
is determined by some external factor. What we call 'envelope' here is called 
'router' or 'rotor- router' by Levine |Lev04j . Holroyd et al. [HLM + 08] . and Levine 
and Peres [LP07| . or 'stacks' by Propp and Wilson [PW98j . or finally 'cards' 
by Diaconis and Fulton [DF91J. The 'commutativity property' is usually called 
'abelian property'. As far as we have been able to find out, the earliest reference 
describing such a type of construction is |DF91j . The general framework consid- 
ered by Eriksson [Eri96j includes both the abelian sandpile model and the card 
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game of Diaconis and Fulton. 

The novelty introduced in this work is twofold. The first thing is that even with the 
sleep envelopes this evolution is still commutative. Finally, and perhaps the core 
of this chapter, the random variable 'amount of visits to each site' increases with 
the initial amount of particles and decreases with the amount of sleep envelopes 
(and therefore it will be probably wise to try and approach the model by looking 
at this variable). 

Before proceeding let us keep at the informal level for one more paragraph and 
explain how these properties may be used. 

First, because of the commutativity property, one has the right to assume the label 
sequence is whatever looks convenient. By doing so one can freely manipulate the 
scenario and dictate the destiny of each particle. The evolution so obtained may 
have nothing to do with real life, but as long as it concerns the number of envelopes 
burned at each site (which is essentially the local time), things are fine. Second, 
because of the monotonicity, the more particles are added to the initial state 
(which happens when the initial box [— M, M] d or the parameter \i increase) and 
the less sleep envelopes are present (which happens when A decreases), the more 
envelopes will be burned. Therefore one can for instance ignore the existence of 
some sleep envelopes when manipulating the scenario, obtaining as a final result 
an upper bound for the local times. All this will be made more clear in the next 
section, when these properties will be used in practice. 

A formal description of our construction is given now. 

Let M be fixed, we shall construct the process having law P^-. For each site 
||x|| ^ M, let t]x(0) be i.i.d. with distribution Poisson(/i). Let m = ^2 x Vx(^) be 
the total number of particles. Label these particles p±, . . . , p m , let xi, . . . , x m be 
their corresponding positions and let T\ — ■ • • — r m = B denote their type, that is, 
'active' or l B\ Take a Poisson process (0 < t\ < t 2 < ■ ■ ■ ) on R + with intensity 
m(l + A). Now take a sequence of independent, uniformly distributed indexes 
N = (ni,ri2, ns, . . . ) where rii G {1, . . . , m}. For each site x G Z d , take an i.i.d. 
sequence (F x j)j eN , where F x j is —1 with probability +1 with probability 
and with probability Let (l be the space consisting on elements u> = 
((x,t),F,T,N) such that x G S m , r G {A, B} m , F G {-l,0,l} 5xN , T = (0 = 
t < t± < t 2 < • ■ • ) with ti t oo, and iV G Af m for some m ^ 0, where Af m C 
{l,...,m} N denotes the set of sequences where each symbol appears infinitely 
many times. Let denote the probability distribution on fl given by the above 
description. 

We now show how to construct the process r]f' B (x), x G S, t ^ from Co. The 
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system state is represented by X = ((x, r), J), where (x,t) = (x n ,T n )™ =1 , J = 
(jx)xes, x n gives the position of the particle p n , r n gives its state (A = passive, 
B = active) and j x tells how many envelopes have been burned at site x. We 
denote the coordinates of X by x n (X), r n (X) and j x (X). The set of possible 
states with m particles is denoted X m = S m x {A, B} m x Z+. 

To recover the original representation, where the particles are undistinguishable, 
we define rj [X] by 

m 

V A (x) = J2^{x n (X))5 A {r n (X)) 

71=1 

and analogously for r] B (x). 

With F fixed, we take Xo = [(x, r), Jo) , where J = and lj = ((re, r), F, T, JV) . 
We shall obtain by updating according to the action performed when the 
clock rings for particle p ni , and we denote this operation by Xi = rti ■ X^\. This 
yields a sequence of states (X , X 1: X 2 , . . . ). 

The process (?7 t (u;)) t>0 will then be given by r)t(u>) = rj[Xi] for t G [U,ti + i). It is 
easy to see that that the process (r)t)t^o so obtained is a cadlag strong Markov 
process with countably many possible states and its infinitesimal generator indeed 
corresponds to the dynamics we are considering, that is, V^(r](uj) E A) = ¥^(A). 

It remains to define n • X. When the clock rings for some active particle it will 
look for the next envelope in the sequence of the corresponding site x, that is, 
F x j + i and perform its instruction: if F x j + i = ±1 the particle jumps to site x ± 1, 
if F x j + i = the particle tries to fall asleep (and succeeds iff it is alone at site x). 

Let n E {1, . . • , m} and X e X m . We are going to define n ■ X. If r n = A, take 
X' = X. Otherwise, let x = x n (X) and j = j x + I. If F- -- = take x' n = x n for 
all n and 

r' n = ^ A, n = h and xi 7^ x n V I 7^ n 

n = h and xi = x n for some / 7^ n 



If Fa j = ±1, take 



and 



x n , n^n 
x n ± 1 , n = h 

Xn 7^ ^ 
■^1 -t 1 
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Take 




3xi x ^ x 

jx 1) 



Finally take X' = {(x' n , r^) n , and put n • X — X'. 

We now state and prove the properties already mentioned. 

It is only assumed that the evolution stabilizes. By that we mean that there 
exist to and r\ such that r] t (Co) = rj V t ^ to, so after some time all particles will 
be passive. We denote by Cl s C Cl the set of Co that stabilize. Though one will 
immediately notice that V^(fl s ) = 1, the argument hold in great generality and 
it involves no probability. 

For Co G f2 s , let rjoo(Co) denote its final state and let R x {Cj) denote the total number 
of jump envelopes burned at the site x during the evolution. 

Theorem 1.2 (commutativity). Let Co = ((x,r), F,T, N) G Q be a realization of 
the graphical construction that stabilizes. Suppose uj' = ((x',r'), F',T',N') G n 
has the same envelopes and initial state of Co, that is, (x', r') = (x, r) and F' — F. 

Then Co' also stabilizes at the same state i q 00 {Cb l ) = rj^LCo). Moreover, the total 
number of particles that visit each site is preserved: R x (Co') = R x (Co). 

The above property, besides being very useful in practice, says that only (x, r) 
and F matter to decide whether Co is in Cl s and to determine rj^ and R x . The 
theorem below states that R x depends on (x, r) and F in a monotone way. To 
be more precise, we write Co' ^ Co when: all particles present in Co' are present 
in Co at the same position; all particles that are active in Co' are also active in Co; 
the envelopes F' can be obtained from F by inserting some sleep envelopes in the 
sequence (F x j)j £N for some sites x. 

Theorem 1.3 (monotonicity). Let Co = \(x,t), F,T, N) G Cl be a realization of 
the graphical construction that stabilizes. Suppose Co' ^ Co, that is, Co' has less 
particles than Co (and less active particles), as well as more sleep envelopes. 

Then Co' also stabilizes and the total number of particles that visit each site does 
not increase: R x (Co') ^ R x (Co). 

The remainder of this section is devoted to the proof of Theorems ll.2l and ll.3l and 
for another lemma that will be used in the sequel. 

Proof of Theorem ll.2L Of course we should only consider the sequence of states 
of the system, regardless of the exact time they happen in the continuous-time 
setting. 
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Suppose we start at state Xq and the first executed action is particle p n , standing 
at site x and performing the action given by F = F x j. It means r n = B, x = x n 
and j = j x + 1. If F is a jump envelope the particle performs that jump, and 
if it is a sleep envelope the particle becomes passive iff x^ ^ i„ Vn 7^ n, but in 
any case the envelope F = F x j is burned. This gives a new state Xi — C • X , 
Co = C(n,x,j), where X\ differs from Xo by the possible change in the position 
or state of particle p n (and the state of other particles that may be activated by 
that change) and by the change j x 1— > j x + 1. 

So, for C(n,x,j) ■ X to make sense we must have x n (X) = x,j = j x (X) + 
l,r n (X) = B. In this case we say that C = C(n,x,j) and X are compatible and 
operator C is called an elementary operator. 

By C we denote a (possibly finite) sequence of operators C = (Ci, C 2 , C3, . . . ), 
which are either elementary operators or the identity operator. We say that C and 
X are compatible if Ck and [(Ck-i o • • • o C\) ■ X] are compatible for all k e N. 
We denote by X k _i the term in brackets and define C • X — (X ,X 1 ,X 2 , . . . ). 
(When X k = X' for all large k we also denote X' by C ■ X and when C is finite 
we also denote by C ■ X the last state in this sequence.) 

For C = (Ci, C 2 , . . . ) and X to be compatible it is necessary and sufficient that 
the following conditions are fulfilled. 

For a fixed n, write for the subsequence of C made of the elementary operators of 
the form C(h,x,j), that is, the operators that correspond to an action performed 
by particle p h . Write C h = (C k(1 , h) , C k ( 2 ,h), ■ ■ ■ ), where C k ^ n) = C(n,x hn ,j hn ). 
The sequences C n , n = 1, . . . , m must satisfy the two properties below. For I > 1, 
write x = x k( i_ hn) , j = j k (i-i, n ), x' = x k (i, n ), f = 3k(i,n)- First, C n describes a 
trajectory, that is, Xfc(i )Tl ) = x n (X) and, for I > 1, one has x' — x if F XjJ - = or 
x' = rr ± 1 if F x j — ±1. Second, if / > 1 and F Xj _,- = 0, that is, particle n has 
previously burned a sleep envelope, or if I — 1 and r n {X) = A, then there must 
be some other particle that will activate it before its next action. More precisely, 
if F x j = 0, there must be either an n' such that x n i(X k n ljn \) = x or k', n', j" and 
x" = x^fl such that k(l — l,n) < k! < k(l,n), C k > = C(n',x",j") and F x i>j» = ±1. 

For a fixed x write C x for the subsequence of C made of operators of the form 
C(n,x,j), that is, the operators that correspond to an action performed at site 
x. Write C x = (C fe(M) , C k ( 2 ,x), ■ ■ ■ ), where C k ^ x) = C(n hx ,x,j ljX ). Then the 
sequences must satisfy x) = j x (X) + 1 and, for I > 1, x) = j(l — 1, x) + 1. 

So, C and X are compatible if and only if all the conditions above are satisfied for 
all n and all x. In this case the final position x n of each particle p n is determined 
by the sequence C n and its state r n is A if and only if the last operator in C n 
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corresponds to a sleep envelope and after this last operator no other particle is 
found in site x = x n . The final value of j x is determined by the sequence C x and it 
is given by j + 1 where the last operator in the sequence C x is C(n, x,j) for some 
n. 

For n G {1, ...,m}, let C = C(X,n) denote the action performed by particle 
n if the system state is X. Namely, C is the identity if r n (X) = A and C = 
C(n,x,j + 1) if r n (X) —B,x n (X) —x,j x (X) — j. Notice that the operation n-X 
introduced just before Theorem II .21 is given by C(X, n)-X. For X = (m, n 2 , . • . ) G 
{1, . . . , m} N and X = X a state, define C = C{X, N) = (d, C 2 , . . . ) by C x = 
C(X , ni), X\ — C\ • X , C 2 = C(Xi, n 2 ) and so on. Write X • X for the sequence 
(Xq, Xi, . . .). Of course in this case C(X, N) and X are compatible. 

Let ii be a permutation of {1, ... , m}. Denote tt(X) = ((x n -i( n ), ^n- x {n))^=i-, (jx)xEs) 
Then vr(n • X) = 7r(n) • vr(X). Notice that r][X] = rj[X] iff X = tt(X) for some vr. 

For N = (ni,n 2 ,n 3 , . . . ) define ir(N) = {n{ni), 7r(n 2 ), 7r(n 3 ), . . . ) , 7r(C) = 
C(ir(n),x,j) for C = C(n,x,j) and 7r(C) = C when C is the identity, 
tt(C) = (7r(Ci),7r(C 2 ),7r(C73),...) when C = (d, C 2 , C 3 , . . .), finally tt(X) = 
(7r(X ), 7r(Xi), . . . ) for X = (X , X 1; . . . ). It is immediate that ir(C) and 7r(X) 
are compatible whenever C and X are so and 7r(C) • 7r(X) = 7r(C • X). Also 
C(tt(X), tt(X)) = tt(C(X, X)) and vr(X • X) = ir(N) ■ vr(X). 

Let 9 denote a shift, that is, 6X = (n 2 , n 3 , n 4 , . . . ), = (X 1; X 2 , X 3 , . . . ), 
and 9C = (C 2 ,C 3 ,C 4 ,...); for X = (m, n 2) n 3 , . . . ), X = (X , X 1; X 2 , . . . ), 
C = {C x , C 2 , C 3 , . . . ). Then (9C) ■ (d ■ X) = 9{C ■ X), (ON) ■ {n x ■ X) = 9(N ■ X). 

We say that the sequence X = (Xo,Xi,X 2 , . . . ) stabilizes if there is X and k 
such that Xfc = X for all k > ko and r n (X) = A for all n = 1, . . . , m. Of course if 
X stabilizes at X then 6X stabilizes at X and 7r(X) stabilizes at vr(X). 

The theorem will follow from 

Claim 1.1. Lei X and (F x j) be given. Suppose the sequence N ■ X stabilizes at 
X for some N G {1, . . . , m} N . Then, for any N' G M m there is a permutation tt 
such that the sequence N' ■ X will stabilize at 7r(X). 

To prove the claim, let us write d(X',X") = Y. X 3^ X ") ~ Jx(X') when X" is 
obtained from X' by applying elementary operators. In this case d(X', X") will 
be the number of such operators. We shall prove the statement of the claim by 
induction on d(X, X). For d(X, X) = we have X = X, so at state X all particles 
are passive and for any sequence X' we have X' • X = (X, X, X, . . . ) . 

Now suppose d(X, X) > 0. We can assume C = C(X, X) = (Ci, C 2 , . . . ) contains 
no identity operators by suppressing the elements in the sequence X that would 
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correspond to such operators. Take N' = {n' l: n' 2) . . . ) G Af m . Since d(X, X) > 0, 
there is at least one n G {1, . . . , m} such that r n (X) = B. Take i 1 = min{i : 
T n >(X) = B} G N, such %\ exists by definition of Af m and the last observation 
(ii is the first time in the sequence N' when the clock will ring for a particle 
that is active in X). Writing C = C(X, N') = (C[, C' 2 , . . . ), we get C[ = C' 2 = 
■■■ = C' il _ 1 = Id and C' ix = C(n,x,j + 1), where n = n' iv x = x n (X), and 
j = jx(X). We can assume %\ — 1 by considering 6 H ~ 1 N' instead of N'. Take 
i* = min{i : Cj = C(n,x,j + 1) for some n,j}, such i* G N exists since C ■ X 
stabilizes, r n = B, and x n (X) = x. Write C { * = C(n*,x*,j* + 1). Since C and X 
are compatible, we must have j* = jx(X) = j because of the definition of i*. 

There are two possibilities. 

First case: n = n*. As we shall see below, if one takes N" = 
(rii*, rii, ri2, . . . , Hi* -i, rii* + i, 7ij*+2, . . . ) then it is the case that N" ■ X stabilizes 
and its final state is given by X. Now n'[ = n[, so N" ■ X = (X,0N" ■ X'{) 
and N' • X — {X,6N' ■ X[), where X[ = n[ ■ X = n'[ ■ X = X'{. But 
d(X 1 ,X) = d(X,X) — 1, so we apply the induction hypothesis for X[, ON" and 
ON' to conclude that N' -X stabilizes at state n(X) for some 7r, finishing the proof 
for this first case. 

Second case: n ^ n* . Take n = ir n ' n * given by 

{n, if n — n* 
n*, if n = n 
n, otherwise, 

and take N'" = (ni, n 2 , . . . , 7r(nj*), 7r(rii* + i), 7r(ni» +2 ), • • • ) • Write N'" ■ X = 
{X%,X'{',X!j!,...). Of course X'g = X k for < k < i*. By definition of i*, 
since N and X are compatible, it follows that x n *(Xi*-i) = x s (Xj*_i) = x and 
Tn.iXr.J = TniX^^) = B, so tt^-.x) = X^. Since O^^N'" = n^^N), 
i *~ 1 X= (^*- 1 X)-X i *_i and i *- 1 X" / = (0 < *- 1 iV" / )-X i ._i = 7r((^- 1 X) 
we have that X'" stabilizes at 7r(X). Also, d(X, 7r(X)) = + 
d^.-ij^X)) = tZ(X, +rf(X i *_ 1 ,A > ) = rf(X,X). So it suffices to show 
that there is some n' such that TV 7 • X stabilizes at n'(X). When one compares 
N' ■ X and N'" ■ X, and repeats the current proof with N'" instead of N, it will 
happen that the first case will hold and, as it has already been shown, there will 
be some n" such that N' ■ X stabilizes at ir"(ir(X)). Take it' = it" on to complete 
the proof of the second case. 

To finish the proof we must show that N" ■ X stabilizes at X in the first case. 
Since e N" = l * N we have l * (N ■ X) = (0 l * N) ■ JQ* and e (N" ■ X) = 
(O^N") ■ X% = (0 e N) ■ X£, so it suffices to show that X£ = X*.. By defini- 
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tion of n* and i*, we have that x k 7^ x* and n k 7^ n* for k = 1, . . . , i* — 1, where 
C fc = C(n k ,x k ,j k ). Take iV = (n x , n 2 , . . . , ?v) and X" = (n", n^', • ■ ■ , n-*) = 
(nj», ni, n 2 , . . . , nj»_i). Let us show that C(X,N") = (Ci*, C±, C2, . . . , Cj*_i). 
First, = n* and by the definition of n* we have x n *(X) = a;*, j x *(X) = j* 
so C(X, n") = C(n*, x*,j* + 1) = Cj* = C". Now X and X" coincide except that 
the position and state of particle p n * may be different, j x * has increased by one 
unit and some other particles in X" might be in state B instead of A because of a 
jump of p n *. So Ci = C(X,nx) = C(ni,x 1 = x ni (X),j 1 = j Xl {X)) = C(X",nx), 
since r ni (X) = T m (X") — B, n\ 7^ n* and x\ 7^ x*. Again, X 1 = C% ■ X and 
X% = d ■ X'{ coincide except that j x *(X'{) = j x *(X) + 1, r n (Xi) = A, and 
TniXip) = B for the particles p n that were activated by the first jump of p n *. By 
repeating this procedure, we get C(X' k ',nk) = C(X k _i,n k ) for k = 1, . . . ,i* — 1; 
X"* and Xi»-i coincide except that j x *(X",) = j x *(Xi*-i) + 1, r n (Xj»_i) = A, and 
possibly r n (X") = B for the particles p n that were activated by the first jump 
of p n *. Now C(Xj*_i,n*) = C(n*, sc*, j* + 1) = Ci*, since a: n *(Xj*_i) = x* and 
j x *(Xi*^i) = j*. So Xi* = Ci* ■ and Xj» differs from Xj»_i on the following 

features. First, an increment at j n * and the change in the position or state of p* n 
according to J x *j* + \. Finally, they may differ on the state of particles that would 
be activated after this last jump, which are all those activated atlw X" except 
those that have been activated and left to some other site during the course of 
X 1— > Xi 1— > . . . h- ► Xj*_i. This implies Xj* = X"„. This finishes the proof of 
Claim 11.11 and of the theorem. □ 

Proof of Theorem 11.31 We have to prove the statement of the theorem when 
Co' G Vt is obtained from Co G Qs by modifying the following features: some 
particles may be removed at the initial state, particles that are kept may start 
passive instead of active, and sleep envelopes may be inserted. We can perform 
such modifications one by one, so it suffices to prove the statement when either 
only one particle is removed, or one particle starts passive instead of active, or 
else one sleep envelope is inserted at a given site. 

We continue to use the objects, notations and results of the proof of Theorem II .21 

Given a sequence C compatible with a state X there is (at least) one sequence 
N such that C = C(X,N). We shall denote by N(C) the smallest one. For 
C = (Id, . . . ,Id, C(ni, x\,ji), . . . , Id, . . . , C(n 2 , £2, J2), ■ ■ ■ )> remove the identity op- 
erators, obtaining C = (C(ni,xi,ji),C(n 2 ,X2,j2),C(n 3 ,x 3 ,j 3 ) y ...). Of course C 
and C are compatible with the same states X and furthermore C ■ X and C ■ X are 
the same except that some positions of the later sequence may appear repeated 
several times on the former one, in particular C ■ X stabilizes at X if and only 
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if C ■ X does. We define N(C) = (n\, n 2 , TI3, . . . ) in this case and it follows that 
N(C) ■ X = C ■ X. 

Let us consider the first case where the envelope configuration is changed so that 
(F X A is given by inserting one sleep envelope at one given site in (F x> j). Take 
N G M m and consider the sequence N - F iX. If the sequence stabilizes before even 
burning this extra envelope, the result follows trivially. So suppose otherwise there 
is some k, n such that C(X k _ 1 , n; F') = C(n, x,j + 1), where FL j +1 is the inserted 
envelope. Then X k is given by X k -i except that the counter jx, increases by one 
unit and that maybe r n (Xk) = A. The first difference does not violate the state 
of the theorem since, for its purpose, burned sleep envelopes are not counted. The 
other difference, that the particle is passive in the evolution ruled by F' when it 
should be active in the evolution ruled by F, is exactly the second case discussed 
below. 

Second (X) = x n (X') for all n, r n (X) = r n (X') for all n = 1, . . . , m except 

for the last particle p m , for which r m (X) = B, r m (X') = A (we have assumed for 
simplicity that the difference is for particle p m ). We have to prove that for some 
N G M m , N ■ X' stabilizes at X' and j x (X') ^ j x (X) V x, under the hypothesis 
that JV • X stabilizes at X. So fix N G M m . Take C = C(X, N) and N = N(C). 
Then N = [n\, . . . ,n^), where k = d(X,X). 

Take N' = (1,2, . . . ,m — 1) and XL-. = N' ■ X', where by abuse of notation we 
call N' ■ X' the last term in the m-tuple. If X',y = X' it means T n (X') = A V n = 
1, . . . , m — 1 and, since r m (X') = A, we have by Theorem II .21 that X' = X' . Then 
of course j x (X) ^ j x (X) = j x (X') = j x (X') and we are done. If T m (X' {1) ) = B, 
since m £ N' we have XL = N' ■ X, so take N' = (N', N) G M m , it follows by 
Theorem 1 1 . 2 1 1 hat N'-X and N'-X' coincide after the first m coordinates, therefore 
N' ■ X' stabilizes at X and we are done. So suppose XL 7^ X' and T m (XL) = A. 
Then d x = d(X',X' {1) ) ^ 1. Take X' {2) = N' ■ X[ x) . If X' {2) = X' (l) it means 
T n (XL) = A\f n = 1, . . . ,m — l and, since T m (XL) = A, we have by Theorem 11.21 
that X' {1) = X'. Then of course Jx (X) > j x (N' ■ X) = j x (N' ■ X') = j x (X>) 
and we are done. If r TO (X( 2 )) / = B, since m G" N' we have XL = N'-X, 
so take N' = (N 1 , N',N) G M m , it follows by Theorem O that N' ■ X and 
N' ■ X' coincide after the first 2m — 1 coordinates, therefore N' ■ X' stabilizes at 
X and we are done. Now if we suppose XL 7^ XL we have d(XL,XL ) ) ^ 1, 
so d 2 = d(X',XL) ^ 2. Now either we will continue to define XL,XL, . . . 
indefinitely, or this construction will halt at some X^)- All we need to do is to 
rule out the first possibility. For that purpose we are going to prove that d k ^ k; 
this suffices, since by the above construction we have d k ^ k. For simplicity we 
consider k = 2. Let C 2 = C(X' , (N' , N')) , N 2 = N(C 2 ) and C 2 = C(X',N 2 ). 
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Briefly, C 2 is just the <i2-tuple of elementary (non-trivial) operators that take X' 
to X'L\ ■ Since the only difference between X and X' is that at state X the particle 
p m is active, C 2 being compatible with X' implies it is compatible with X and 
C 2 ■ X = N 2 ■ X. By Theorem OJ (N 2 , N) ■ X stabilizes at X, so d(X, X) ^ d 2 , 
finishing the proof of the second case. 

Third case: X = ((x n ,r n )™ =1 , (j x ) xeS ), and X' = ({x n ,T n )™=l, (jx)xes) (we have 
assumed for simplicity that the missing particle is p m ). As before, we have to 
prove that for some N G M m -ii N ■ X' stabilizes at X' and j x (X') ^ j x (X) V x, 
under the hypothesis that N • X stabilizes at X. We fix iV G Af m and take 
C = C(X,N) and N = N(C). Then N = (m, . . . , n k ) , where k = d(X,X). We 
consider the same N', and keep on defining XLyXL^,... as before, and now 
the condition for this construction to stop is when X', k -, = X'^ k _ 1 y which means 
T n (X', k _y) — A for n = 1, . . . , m — 1. Again, if Ck is compatible with X 1 , it is also 
compatible with X, with some abuse of notation since now we are considering 
elementary operators acting on different state spaces. So by the same argument 
as for the second case, k ^ dk = d(X',X', k -,) ^ k and for some finite k we have 
X' = X[ k) = (C' dk o...o C[) ■ X' satisfies r n (X {k) ) = A for n = 1, . . . , m - 1. 
Again in this case X^) = (C' dk o ■ • • o C[) • X is well defined and it differs from 
X' for the presence of an extra particle and because for some n it might be that 
r n (X {k) ) = B and r n (X') = A but j x (X')= j x (X (k) ). Now j x (X (k) ) ^ j x (X) since, 
by Theorem ll.2[ • Xr^) stabilizes at X. This completes the proof of the third 
case. □ 

A different coupling allows the comparison between different initial conditions 
for a fixed time. We shall need the fact below for a mild technical question of 
exchanging the order of a certain limit in Section II .41 For fixed time t > let R x 
denote the number of jump envelopes burned at site x up to time t. 

Lemma 1.1. Let P,P' denote probability distributions on E" such that P >z P' . 
Take V and V as the corresponding laws on Q according to the construction 
described above. For each x G S, t > 0, and r G N, V(R X ^ r) ^ V'(R X ^ r) . 

We shall just highlight the main steps of the proof: 

Let (j)o' B (x)) and (F x j) x j be given. We shall construct (r) t )t^o i n a slightly 
different way. Take (T x j) x j as i.i.d. exp(l + A) random variables and T x j = 
J2{T x ,i- Let J° = for x G S. Write D x {t) = f*rjf(x)ds. Take t x = mi{t : 
D x (t) = T x j = i(t) for some x G S} and take X\ as the (a.s. unique) x such that 
Dx(t\) = T X j = i{t\). Let rj s = i] Q for s < t\ and i] tl be given by the new state 
obtained from r] by opening envelope F xis i and set J], = S xi (x). 

Suppose rj s has been defined for s G [0, t k ]- Take t k +i = mf{t : D x {t) = 
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T x J k + i{t) for some x G S} and take x k +i as the (a.s. unique) x such that 
D x {t k+l ) = T x>J k +l {t k+1 ). Let 7] s = 7] tk for t k ^ s < t k+1 and r] tk+1 be given 
by the new state obtained from r] tk by opening the envelope F x tJ k +1 . Set 
J x + = J x + 5 Xk+1 (x). 

Now by a routine argument, similar to that used to prove super-additivity for first 
passage percolation, one proves the stochastic domination claimed on the lemma. 
□ 

1.4 The phase transition 

In this section we present the proof of Theorem ll.lt The proof strongly relies on 
the commutativity and monotonicity properties shown in Section 11.31 and approx- 
imations by finite boxes (11.21) will be important. 

We shall first give a simple characterization of phase transition in terms of burned 
envelopes. From this we immediately prove the existence of \x c G [0, oo] and that 
fi c is increasing in A. Then we show that fi c $C 1 for one-dimensional simple 
random walks. Finally, in the remainder of the section we give the proof that 
/i c ^ for simple random walks, which is the main result. 

By translation invariance, the phenomenon of fixation is equivalent to the fact 
that a.s. there is fixation at the origin, that is, for large enough time % (0) = 0. 
Consider the random variable R\ given by the number of times a particle jumps out 
of the origin up to time t. As t increases one has Rq f Rq and fixation is equivalent 
to R Q < oo. Define the event A\. = [Rq ^ r]. Now A\ tt^oo A r | r ^oo A^, where 
A r = [R Q ^ r], and fixation is equivalent to P M (A 00 ) = 0. 

Note that by ( II. 2ft one has 

W{Al) = lim F» M (Al), 

M-+oo 

so fixation means 

P^(A r ) — ► uniformly in M. (1.3) 

r— >oo 

Now, within the framework of Section 11.31 the event A r means that at least r 
jump envelopes are burned at the origin. 

It is immediate that (II .3p is sufficient for fixation. To see why it is necessary we 
remark that P^ f (A*) is increasing in both M and t, so switching limits we get 
P^(A*) ^ P M (A r ). Monotonicity in t is obvious and monotonicity in M follows 
from Lemma [1. II 
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Now suppose fl X . 3 j) holds for some /i. By Theorem II .31 the same limit will hold for 
any // < /i. So one can take \i c = sup{/z ^ : limit (11.31) holds}. (Notice that it 
does not exclude the possibilities /i c = or /i c = oo.) It follows from Theorem 11.31 
that /i c is nondecreasing in A. All this holds in any dimension d ^ 1. 

To prove \i c < oo in one dimension we argue as follows. Suppose fi > 1. With 
high probability there will be an amount of particles of order fj,M in a given large 
region (say, in [— M, 0]). Since no more than M particles can stay passive in this 
region, an excess of particles (at least about (/i — 1)M) will have to exit such 
region, and because of the topological constraints of the one- dimensional lattice, 
they must exit either by its left border or by its right border, that is, the origin. 
This implies that whp at least (/i — l)M/2 particles will eventually pass by the 
origin and therefore Rq = oo almost surely, so the system does not fixate. 

This actually proves that fi c ^ 1. We have implicitly assumed that the random 
walks are of nearest-neighbor type but with slight modifications in the argument 
one easily shows that /i c ^ 1 for any finite-range random walk. 

Let us now show that [i c ^ ^ or nearest-neighbor walks. With simple adapta- 
tions the same argument proves that /x c > for finite-range random walks. 

First, we give a heuristic proof. The rigorous argument follows next. 

We shall build a global trap for the particles using disjoint sites to guarantee that 
they cannot collaborate (activate each other) to overcome the trap. The density 
of individual traps one is going to find per site will be given by because this 
is the probability that between two consecutive jump envelopes there is at least 
one sleep envelope. If the density of particles is smaller, that is, \i < then on 
a large scale the system will be trapped whp. In order to be able take advantage 
of this large-scale behavior, we choose a large enough region A around the origin 
and we sweep A. Once this is done the particles cannot come close to the origin 
again because, due of the large-scale domination of the traps over the particles, 
the probability that the trap fails will be small. 

To 'sweep' the region A we let the particles move inside this area until they reach 
its boundary. The particles will be swept one by one, this is possible because of 
Theorem 11.21 None of them will become passive until this step is finished (here 
we use Theorem 11.3!) . During this procedure the probability that more than r 
particles cross the origin is small if r is taken large enough. 

We prove now that whp no particle will come back to the origin after A is swept. 
Place an imaginary barrier at the origin and pick the particle that is closest to 
this barrier. Let this particle jump (Theorem II . 21 again) until it hits the origin (or 
some fixed site far away). While we make the particle jump we are ignoring that 
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each envelope this particle is opening could have been a sleep envelope (Theo- 
rem H3] again) . Then we start coming back along this path while looking for sleep 
envelopes that might have been missed, until we succeed to find the first one. 

The set of envelopes we have discovered so far has the following property. Suppose 
the label sequence were exactly as we had assumed, and the sleep envelopes we 
have ignored were indeed absent. Then all the particles that started in A would 
move until they reach the border of A. Next, one of these particles would start 
jumping until it opens the sleep envelope we have managed to find. Therefore this 
particle would become passive before reaching the barrier. 

Now place a new imaginary barrier at the site where the trap for the first particle 
has been set. Pick the particle that is closest to this new barrier. Again let it 
jump until it hits the new barrier. Then trace back its path looking for sleep 
envelopes that had been ignored, until the first one is found. As in the previous 
paragraph, the evolution we are building makes the particles evacuate A, attracts 
the first particle to the first trap and now it attracts this second particle to this 
other trap without touching the origin and without activating the first particle. 

Again place a new barrier where the second particle should be trapped and keep 
repeating this procedure until a trap has been set for each of the particles. The 
global configuration so obtained is enough to guarantee that all the particles will 
be passive without coming back to the origin. 

As mentioned in the previous section, this evolution may have nothing to do with 
that governed by P^. Nevertheless the real evolution (with a random sequence 
of labels and a few more sleep envelopes) will have the same value for R by 
Theorem 11.21 or maybe even smaller by Theorem 11.31 

It is possible that this procedure fails. It will fail if, at some stage, say the n-th 
stage, we do not manage to find a sleep envelope between the position of the 
n-th particle (after sweeping) and the (n— l)-th barrier. But the probability that 
this will happen is very small because the density of traps is smaller than that of 
particles and A has been chosen large enough. 

The remainder of this section is devoted to formalizing the algorithm just de- 
scribed. 

For each x e S split the envelope sequence (F x j)j into two sequences: one con- 
sisting on the jump envelopes and one consisting on the sleep envelopes. The 
first sequence is given by F X) k = where j(x,0) = 0, j(x,k) = min{j > 

j(x, k — 1) : F X) j = ±1}. So F xk is the k-th jump envelope at x and between the 
(k — l)-th and k-th jump envelopes there are D xk = \j(x, k) — j(x, k — 1) — 1] 
sleep envelopes in (F Xs j)j. Of course (F Xt j)j is in bijection with (F Xj k, D x ,k)k- Also, 
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inserting a sleep envelope in (F x j)j means increasing the value of D Xjk for some 
k. 

It is routine to check that (F Xjk ) Xjk are i-i-d. having value +1 with probability p 
and —1 with probability q = I— p. Also (D X)k ) Xyk are i.i.d. geometric variables that 
are positive with probability We finally remark that (F x>k ) x>k and {D Xtk ) Xtk 
are independent. 

We are going to show that the system fixates if p, < j^. We first present the 
proof for the symmetric case p — q — \ and after we give the modifications in the 
case p—\ — q>\. 

Let e > 0. Take K with the following property: if S n and S' n denote the sum 
of the n first random variables of an i.i.d. sequence distributed respectively as 
Bernoulli(Y^) and Poisson(/x)) then the probability that S' n+1 > S n for some 
n ^ K is less than e. The reason for this choice of K will be clear later. 

In the sequel we shall consider the evolution given by the envelope sequence (F x j ) 
instead of (F x j), where (F x j) * s obtained from the same (F x j) but (D' k ) = 0. 
It means that we take the configuration u>' = ((x n ,T n )n,F',T,N') y u>. Here F' 
is obtained by removing all the sleep envelopes and N' will be constructed in an 
algorithmic way. Of course for this strange configuration uj' the particles will never 
become passive since all sleep envelopes have been removed. So later on we shall 
put back a few sleep envelopes present in uj to construct an intermediate scenario 
uj" y uj such that the particles become appropriately trapped in the evolution 
given by uj" . 

Let mf be the number of particles in uj and relabel them p 1 , . . . , p m , . . . , p m > in 
a way that \x n \ ^ K for n = 1, . . . , m and \x n \ > K for n = m + 1, . . . , m' . 
Let X = ((x n , T n ) n , (j x )x) be the initial state, with j x = 0. Take n\ — 1 and 
X\ = rii ■ Xq. If xi(Xi) = ±K we say particle p\ has been swept and we move 
to particle p 2 ; otherwise we take n 2 = 1, X 2 = n 2 ■ X 1 and so on. Once we have 
xi(X kl ) = ±K we move to particle p 2 . Then we take n kl +i = 2, X kl +i = n^+i-X^ 
and keep doing this, thus obtaining X , Xi, X 2 , . . . until x 2 (X k2 ) = ±K for some 
k 2 . After this we take n k2 +i = 3, A fc2+1 = n fc2+1 • X k . 2 and keep doing this 
until X3(X k:i ) = ±K for some k%. At this point particles pi, p 2 and p% have 
been swept to the boundary of [-K, K\. We repeat this procedure until we have 
swept all of pi,..., p m , thus getting (m, n 2 , n 3 , . . . , n kl ,n kl+1 , ...,n k2 ,..., n km ) 
and (X ,X 1 ,...,X km ). 

After this step the region [-K, K] is clean and \x n {X km ) \ ^ K for all n = 1, . . . , m'. 
Moreover, the number of particles at each boundary i] B (±K) with r] = rj[X km } is 
Poisson(/i(fr + 1)) and is independent of (r] B (x))\ x \ >K - Also, j x (X km ) is a.s. finite, 
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so for there is r such that the probability that j x (Xk m ) ^ r at the origin x = 
is less than e. 

We are now going to build a trap to show that after sweeping no particle can come 
back to the origin. 

To start this second step relabel the particles in X^ m , that is, take X = 7i(Xk m ) 
for some ir such that K ^ x\ ^ x 2 ^ ■ • • ^ x m ii and — K ^ x m "+i ^ Xm"+2 ^ 
• • • ^ x m > in X . 

Let a = 0, 6 — [(3 + 1/A)M] and W = 0. We shall build a sequence of intervals 
[a n , b n ] D [a n -i, 6 n -i] 5 • • • where the particles p n , . . . , p m » will remain confined. 
Finally W C C • • • will be the sequence of traps where each particle will 
become passive in a modified configuration Co" -< Co. 

Take n[ — 1, X± — n[ ■ X , n' 2 — 1, X 2 — n' 2 ■ X±, n' 3 — 1, X 3 — n' 3 ■ X 2 , and so on 
until for some l± > the position of particle p\ is either a or 6 - 

The idea is to set up a trap for p\ strictly between the origin and pi's initial 
position in a way that (i) this particle does not reach the origin and (ii) this 
operation consumes as little space as possible. The first requirement is exactly 
what we want to prove and the second requirement will be important because 
when setting up the trap for the other particles p 2 , p 3 , ... we want to guarantee 
that none of them will activate pi again. 

Suppose = a . The set T\ = {a + 1, . . . , xi(X ) — 1} contains the possible 

sites where the trap for pi will be set up. For x G T\ let h\{x) = j x (Xi 1 ). Now 
take x\ = min{a; G T\ : D x>kl ^ > 0}. If this set is empty we have failed in finding 
a trap for p 1 and we stop the construction. Otherwise we declare the first trap 
successful. The construction of the first trap is complete on defining a\ = xi, 
h = b and W 1 = W U {{x u k^))}. 

In the case Xi(X h ) = b we take Ti = {M + 1, . . . , b — 1}, Xi = max{i G T x : 
D x ,ki{x) > 0}j a i — a o, b\ = x\ and the same expression for W\. 

Now let us trap p 2 . Again p 2 should neither touch the origin nor activate pi 
and the trap for p 2 should consume as little space as possible. Take n' li+1 = 2, 
Xi 1+ i = n' h+1 -X h , n' h+2 = 2, X h+2 = n' h+2 - X h+ll n' h+3 = 2, X h+3 = n' h+3 -X h+2 , 
and so on until for some l 2 > the position x 2 (Xi 2 ) of particle p 2 is either a x or 

If x 2 (Xi 2 ) = ai we take T 2 = {ai + 1, . . . , x 2 (Xq) — 1}, for x G T 2 we set k 2 (x) = 
j x {Xi 2 ) ) then take x 2 = min{x G T 2 : D x ^ 2 i x \ > 0}. Again if such set is empty 
we have failed in finding a trap for p 2 and we stop the construction. Otherwise 
we declare the second trap successful. Finally we define a 2 = x 2 , b 2 = b\ and 
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W 2 = W l U {(x 2 , k 2 (x 2 ))}. If x 2 (X h ) — bi we take T 2 — {M + 1, . . . ,bi — 1}, 
x 2 = max{x G T 2 : D x ^ 2 ^ > 0}, a 2 — a\, 62^2 and the same expression for W 2 . 

In this way we construct [ao, bo], Wq, [ai, bi], W\, . . . , [a m », b m "], W m ". Starting 
with ao = — bo, bo = and Wq = we perform a similar construction for 
[d , b }, W , 61], Wi, . . . , [ },w m , 

—m" • 

Suppose all the traps have been successful. Take W = W m "UW m >- m ", the total set 
of traps. W has exactly m! elements, one trap corresponding to each particle. By 
construction D X:k > for all (x, k) G W. Take D" x k = H( x ,k)ew so that F)" xk < D x k 
and consider the envelope sequence {F X J) obtained from {F X J) and (D" k ). Write 
n k m +i = 7r _1 «) for % = 1, . . . , l m > and take N" = N' = (m, n 2 , . . . , n km+ i m ,). Let 
(X^XC, ...) represent the evolution of a)" = ((x n , r n ) n , F", T", iV") b a). We 
leave for the reader the task of checking that x n (X") 7^ for all % ^ k m and that 
the evolution (X", X 2 ' , X% , . . . ) stabilizes at tt~ 1 (X), where x n (X) = x n . In other 
words, the particles leave (-K, +K) in the first step and during the second step 
they get stuck at the traps at x n without being at the origin again. 

The procedure may stop at some point, namely if one fails to set a certain trap. 
Given that the n-th trap was successful and it was placed at site x n , the probability 
that the (n + l)-th trap will be at site x n + 1 + k is given by ( j^) k j^- By the 
definition of K, the probability that one of the m" first traps is not successful is 
less than e. Similarly, the probability that one of the m! — m" other traps fails is 
also bounded by e. By the choice of ro we have that P^(A r ) < 3e for any r ^ ro, 
where ro does not depend on M, finishing the proof for the simple symmetric 
random walk. 

For the simple asymmetric random walk the proof has slight adaptations that we 
indicate in the following three paragraphs. 

First we explain how to define a suitable K. Suppose p > q, so the random walks 
are biased to the right. For k > 0, consider the following system of random walks. 
Start with Poisson(2A;/i) number particles at site x = k and, independently of 
that, Poisson(/x) particles at every site x > k. Then let all these particles perform 
independent random walks with probability p of jumps to the right and q = 1 — p 
of jumps to the left. Let Xq be the random variable given by the minimum position 
among all such paths. It is easy to see that xq — > 00 in probability as k — > 00. So 
choose ki such that x is positive with probability at least 1 — e. We are now going 
to choose k 2 . Let S n (resp. S' n ) denote the sum of the n first random variables of 
an i.i.d. sequence distributed as Bernoulli(Y^) (resp. Poisson(/x)). Then there is 
k 2 > such that the probability that S' n+k > S n for some n ^ 1 is less than e if 
k ^ k 2 . Take K = fc x V k 2 . 
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There is also a change in the sweeping procedure defined on p. [2H Instead of 
stopping to move a particle until it reaches ±K we move it until it hits +K. This 
will happen almost surely since the random walk is biased to the right. 

The second step, namely that of finding traps, is the same. Once K is chosen, the 
choice of r is also analogous. So, the probability that this procedure will fail will 
again be bounded by 3e. 

This finishes the proof of Theorem 11.11 □ 



1.5 Concluding remarks 

The proof that /i c ^ y+a f° r simple random walks presented in Section 11.41 can 
be easily adapted for finite-range walks to show that /i c > 0, though the authors 
conjecture that /i c ^ y+a should still hold in this case. 

The initial distribution being Poisson with parameter /i plays no special role in 
the proof. Indeed the proof of the lower bound for /i c works with any family 
{-P/i}o^i<oo of distributions on {0,1,2,...} parametrized by the first moment. 
For the uniqueness of the phase transition we only need the natural assumption 
that P M >z Pfj,' when /i ^ /i', that is, the bigger /i, the more particles. 

Yet we conjecture that this phase transition is actually universal with respect 
to the distribution. There should be a /i c such that for any spatially-ergodic 
initial distribution the system fixates if E[r] AB (0)] < fi c and does not fixate if 
^[^^(O)] ^ ^ _E[?7 A (0)] > 0. In the later case the system should converge to an 
ergodic distribution that depends only on /i. 

Some questions remain open. We do not know how to answer such questions, but 
nevertheless we strongly believe that the same scheme developed here could be 
extremely useful for this purpose. 

To complete the qualitative description of this model, it is still open to rule out two 
possibilities, which are conjectured not to hold. These possibilities are physically 
absurd, but so far there is no mathematical proof that they do not hold. We have 
shown that \i c increases as A increases, that < \i c ^ 1 for all A and that fi c — > 1 
as A — > oo. What remains to be proven is that fi c < 1 for all A and fi c — > as 
A — > 0. Strikingly enough, even the fact that for some A > and /i < 1 the system 
does not fixate (which would be implied by either of these) remains unproven. 

One can also consider /x fixed and study the phase transition at A c , as was exposed 
at the introductory chapter. Theorem 11.11 then says that A c < oo for \i < 1, that 
A c = oo for n ^ 1 and that A c — > as fi — > 0. What remains open then is that 
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A c — > oo as /i y 1 and that A c > for all jj > 0. 

In dimensions d ^ 2, the fact that /i c > is a very interesting open problem. 
Whereas for d — 1 the trick of letting the particles evolve until hitting barriers 
and then setting traps worked well, for higher dimensions different ideas will be 
necessary. 

The proof that \x c ^ 1 for dimensions d ^ 2 was recently given by Shellef [She] . 
his proof is quite simple within the framework we consider here. The fact that 
/i c < 1 still remains as an open question in any dimension. 

The open problem that is, in my opinion, the most interesting but the hardest to 
approach is that of whether there is fixation at \x = \i c . Unlike the most common 
models that exhibit phase transition, such as percolation or Ising model, we believe 
that the critical behavior of this model bears the most important features from 
the supercritical phase. Namely, at // = \i c the system should not fixate. The 
only case considered so far [HSj is that of one- dimensional, totally asymmetric, 
nearest-neighbor random walks. 



Chapter 2 

The Broken Line Process 



In this chapter we introduce the continuum broken line process and derive some 
of its properties. This process is similar in spirit to the Arak-Surgailis process. 
Its discrete version, introduced by Sidoravicius, Surgailis and Vares, is presented 
first. A natural generalization to a continuous object living on the discrete lattice 
is then proposed and studied. The broken lines are related to the Young diagram 
and the Hammersley process and are useful for computing last passage percolation 
values and finding maximal oriented paths. 

For a class of passage time distributions there is a family of boundary conditions 
that make the process stationary and reversible. For such distributions there is a 
law of large numbers and the process can be extended to the infinite lattice. 

One application is a simple proof of the explicit law of large numbers for last 
passage percolation known to hold for exponential and geometric distributions. 

This chapter is based on a research program carried out in collaboration with 
V. Sidoravicius, D. Surgailis, and M. E. Vares. 

2.1 Introduction 

The main motivation of our approach is the discrete geometric broken line pro- 
cess. Informally it can be described as the following particle system with creation 
and annihilation. Consider a sequence £ t ,x of independent and geometrically dis- 
tributed — k) — (1 — A 2 )A 2fc , k = 0,1,... with parameter < A 2 < 1) 
random variables, indexed by points (t, x) e Z 2 such that both t and x simultane- 
ously are either even, or odd (the set of all such (t,x) is denoted by Z 2 .) Assume 
that at each point (t,x) G Z 2 , £ tjX pairs of particles are born with opposite veloc- 
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ities ±1. The born particles then move with constant velocities. At the moment 
when a moving particle collides with another moving particle having opposite ve- 
locity, they annihilate each other. Because of the possibility of multiple particles 
being at the same time-space point, the annihilation rule says that at the moment 
of collision, first annihilated are 'older' particles. 

Trajectories of particles in space-time give rise to a process of (partially coincid- 
ing) 'discrete' random broken lines whose segments (corresponding to individual 
particles) lie at the angle ±7r/4 to the time axis. The broken line process is 
stationary in time and space with respect to natural translations of Z 2 and can 
be defined on Z 2 by consistency, as the limit of similar broken line processes on 
bounded domains. Let <f>t )X be the number of broken lines which pass through 
the point (t, x) G Z 2 . The process {4>t,o} is called the intersection process, it is 
stationary but not i.i.d. and it satisfies a law of large numbers. This setup is 
described in Section |2T2| which is based on the unpublished preprint [SSVj . 

A generalization of this object is what we call the continuum broken line process, 
introduced in Section 12.31 Within this framework, instead of having a certain 
number £ G Z + of pairs of particles being born at each space-time point, we 
have a mass whose value is in principle non-integer. A few subtle differences will 
arise from this difference but one can overcome them by looking for broken lines 
of finite size, even when the process is defined on the whole lattice. In view of 
Proposition 12.11 and its explicit construction, it will hopefully be clear that the 
broken line process and its continuum version are indeed natural objects. Also 
presented in Section l2~3l is the discussion of necessary and sufficient conditions for 
reversibility, which turns out to be crucial for the applications that appear after. 
Applications to last passage percolation are shown in Section 12.41 



2.2 The geometric broken line process 

In this section we shall describe the evolution of a particle system with discrete 
time step. This system is first defined for finite space-time domains. Later on it is 
described as a non-homogeneous Markov chain and we shall see that this evolution 
is reversible and consistent. As a consequence of the above properties, the model 
is extended to the whole lattice. We shall define the broken lines starting from the 
space-time trajectories of the particles. It is instructive to consider the discrete 
geometric broken lines process first, before proceeding to full generality. This 
presentation is based on [SSVj. 
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2.2.1 The evolution of a particle system 

Let 

S = {(t, x) e Z 2 : t ^ t ^ h,x t - ^ x < x t ,+ }, 

where to,ti G Z, t < t\ are given points, and x t - ^ x t,+ ,to ^ t ^ t\ are paths 
in Z 2 such that for some t ^ ^ ti, £t+i,± — x tj ± = ±1 (t ^ t < t^) and 
x t+i,± ~ x t,± = (^i ^ t < ti); in other words, x t)+ increases on the interval 
[ioi^i] an d decreases on while x t _ decreases on [to>^i] an d increases 

on [t ,i)^i]- Any such set S C Z 2 will be called a hexagonal domain. Let 



Figure 2.1: a hexagonal domain 

x to - ^ x ^ x t ,+ } denote the left vertical boundary 
of S, and d+S = {(t,x) G S : t ^ t < tj^x = x ti+ },d_S = {(t,x) 6S:i < 
i ^ ^o,i5 x = ^t,-} be the 'northwest' and the 'southwest' boundaries, respectively. 
Put also d° + S = d + S\d S, d°_S = dS\d S. 

For any point (t,x) G S, let r)£ x ,r)^~ x denote the numbers of ascending (moving 
with velocity +1) and descending (moving with velocity —1) particles which 'leave' 
(t,x), and Ctx) Ct~x the respective numbers of ascending and descending particles 
which 'come' to (t,x). The relation 

Ctjz ~~ Ct,x = Vt^x ~ Vt,x (2-1) 

says that at each point (t, x) G S, particles can be created or killed by pairs only. 
More precisely, we assume that each pair of incoming particles with opposite 
velocities is annihilated, and a certain number £ tjX = 0, 1, . . . of pairs of particles 
with opposite velocities is born at (t,x), so that 

= + [C£ - + , %x = it,x + [Ci, x - Cjj + - I 2 - 2 ) 



32 



CHAPTER 2. THE BROKEN LINE PROCESS 



Moreover, 

C£ = vt w (t,x) e S\(d T S U d S), (2.3) 
as all transformations of particles may occur only at lattice points (t, x) G S. Put 

V = (v£x,V^x ■ (t,x) G S), 

C = (CLC t : x -(t,x)ed S), 
( + = ((; x :(t,x)ed° + Sl 
(„ = (C+ x :(t,x)ed°_S). 

(° = (Co, C+)C-)> an d regard 77 as configuration in S and (° as boundary data. 
It is clear that there is a 1-1 correspondence between (C°>£) and where 
(£ x ,rf x ,£t,x are related by (|21jl - (|23|l . Note that (fill]) and (Q imply 

- = ryti^-i - %x+i, (t, x) G S\(d S U d+S UdS). (2.4) 

The probability measure Ps{v) corresponding to the evolution of this particle 
systems in the hexagonal domain S can now be defined as follows. Let < A < 
1 be a parameter. Assume that all CtV s m ^° are independent and Geom(A)- 
distributed, i.e. 

p(£± = k) = (1 - A)A fc , k = 0, 1, ... , (t, x) G <9 .S U <9 T S. 

Assume also that all £t i!r 's in £ are independent and Geom(A 2 )-distributed and, 
moreover, (° and £ are independent. Then Ps is the distribution of 77 = {rjf x) r)^~ x : 
(t,x) 6 S). In other words, Ps is the distribution of the number of outgoing 
particles from (t, x) G S, assuming that C^ xl (t, x) G 8qS U d T S (the number of 
particles which 'immigrate' to S through its left boundary d S U d + S U d^S) are 
i.i.d. Geom(A)-distributed, and £ tjX ,(t,x) G S 1 (the number of pairs of particles 
with opposite velocities born inside S) are i.i.d. Geom(A 2 )-distributed, indepen- 
dent of the C^a;' s - in the evolution, the born particles move with constant velocities 
+ 1 or —1 until they collide at some lattice point (t,x) G S, after which the col- 
liding particles die (annihilate). 

2.2.2 An equivalent description for the evolution of the 
particles 

Below we provide another description of the evolution of particle systems intro- 
duced above. In this setting various properties are more easily seen to hold true. 
Let Eg be the set of all configurations 77, where r\ t x take values 0, 1, . . . and satisfy 
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relations ( 12.41) . Let £d°s be the set of all configurations (° = (Co, C + ? C— ) ? Co = 

(Ck,Ct7» : e = (CtT* : e ^),C- = (C£ = (M) e d°_S) on 

c^S 1 = <9oS'U<9 + S'UcLS'. Let £5 be the set of all configurations fj = 77), where 
C° G T.Qog,!] G S5, and, moreover, 

{Cm - ^i,x+u if (*» x ) G d-S\d S 
- Cx, if (*> e d + S\d S . (2.5) 
C^-CiTx, if (*, x) G 5 5 

Finally, for a given configuration £° G let Sg^o denote the class of all rj G £5 

which satisfy ( 12.51) . 

We shall define a probability measure Qs\c°(v) 011 ^s|C° as a ( non homogeneous) 
Markov chain r] t ,to ^ t ^ t\, whose values at each time t are restrictions of a 
configuration r\ = {v^ x) rj~ x : (u,x) G S) G Ss|c° on i = u; in other words, 
Vt = (vt,x,Vt,x ■ x E S t ), where S t = {x : (t, x) G S). 

The transition probabilities of the Markov chain are defined as follows (for clarity, 
we assume tgi = ^01 = ^01 below). 

(i) At t — t , the distribution of i] to = (j]f o X) rj^ x : x G doS) depends only on 
Co = (Ci> Qx : (*, x ) e <9 >S) and is given by 

<55ic°(%,ico) = n ^i^'^ici^CtTJ, (2.6) 

where 

g(n + ,n"|m + ,m") = Z~+ m _A n +n 5 {n +_ n - =m +_ m - } , (2.7) 
n*, = 0, 1, . . . , < A < 1 is a parameter, and 

^ ^|m + — m~\ 

(ii) Let to < * ^ £01 • Then = (r)^ x ,T]f X ,x G St) depends only on ^ t _i = 
(vt-i,xiVt-x,xi x G St-0 and Ct,+ = C^x ti+ ,(t - = Cte*,-i according to the transition 
probability 

Qs\c(vth-i,(t,±) = q(v^ t , + ,vr >X t,+\vti,x t - 1 ,+X^ xt>+ ) 

x i(vt,x t ,-,v7,x t ,Jvti,x t -i,-it£x t ,J ■ (2.9) 

X Y[xdSt,x^xt,± Qv)t,xi T )t,x\ T )t-l,x-li Vt-l,x+l) 

(iii) Let i i < t ^ t\. Then r\ t = (T)txjVixj x ^ depends only on 77^—1 = 
i,xiVt-i,x> x ^ St—i) an d 

Qs\t°(vth-i) = IJ ^fe'^l^-M-i'^-M+i)- ( 2 - 10 ) 

x&St 
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Let IIa(C ) be the product geometric distribution, i.e. all random variables 
(£ x , (t,x) G d S U d-S,C£xi (tj x ) £ 9 S U d + S are independent and Geom(A)- 
distributed. Define the probability measures Qsijf) = Qs(C°^ r l)jQs( r l) on ^s, 
by 

Qs(C,v) = MC)Qs\c(v), Qs(v)= E Qs(C,v), 

C°es S o S 

respectively. We claim that 

Qs = Ps, (2.11) 
where Ps was defined in Section 12.2.11 above. Because of the 1 — 1 correspondence 
(C°,v) «-> (0S0> ( 2 - 13 ) follows from 

Qs|C°(0 := Qs|c°(^ = ^C°,v) = £) = n A»(0- (2-12) 
Equation (12.121) can be shown by induction in t\ — t = 0, 1, . . . . Let t\ = t . As 
26,x = Vt + % x - Ki - CJ, (*, e d Q S, see (O), so by flZZ^flZS) 

P(^ = A;|C£: = ^ + 5 C t 7 a> = m_) = 
= g(n + , n~ |m + , m~) 

n++n — =2fc— |m+— m — | 

7— 1 \ ^ \ n+ + n ~A 

~~ ^m+.m- n +_ n - =m +_ m - 

n++n _ =2fc— |m+— m~ | 

= (l-A 2 )A 2fc = vr A2 (A;), 

where 7r A stands for Geom(A)-distribution. Hence and by (12. 6p it follows 
that d2~T2l holds for * = t x . In the general case, (T2~l"2l and (f2~TT|) follow by 
induction in t\ — to an d a similar computation, based on the explicit form of 
transition probabilities Qs\c°(Vt\Vt-i) given by (I2.9p - (l2.10p . 

In the sequel we discuss some properties of Ps(v)- 
2.2.3 Duality 

The transition probability (12. 7p satisfies the following relation 

7r\(m + )7r\(m~)q(n + , n~\m + , m~) = 7T\(n + )iT\(n~)q(m~ , m + \n~ , n + ), (2.13) 

m ± ,n ± = 0, 1, ... , where ii\ is Geom(A)-distribution. A similar duality relation 
holds for transition probabilities of the Markov chain rj t , to ^ t ^ t\. For example, 
if t i < t < ti, then from (I2TT3D . ( I2TT0D it follows that 

)q{j] 

t,xi Vt,x\Vt-l,x-li Vt-l,x+l) 



2.2. TEE GEOMETRIC BROKEN LINE PROCESS 



35 



The above duality implies that the construction of Ps{fj) = Qs(v) can be reversed 
in time. Namely, in the 'dual picture', 7]f x is the number of 'descending' parti- 
cles which 'come' to (t,x) from the right, and (f x is the number of 'ascending' 
particles which 'leave' (t, x) in the same direction. A similar 'dual interpreta- 
tion' can be given to % x ,(tx- m the dual construction, r]° = {{j]f x1 rj~[ x1 (t,x) G 
diS), {rfix-, (t,x) G d + S\d x S), (% x , (t,x) G d~ S\d\S)} is the boundary condition; 
here, d\S is the right vertical boundary, and d + S,d~S are the 'northeast' and 
'southeast' boundaries of S, respectively. The probability measure Qs\r/°(() is de- 
fined on the set Sg^o of all configurations ( = (CtxiCtxi (^> x ) ^ which satisfy 
conditions analogous to (I2.4p . (12.51) ; namely, 

C - Ci = cr + i,,-i - (*, x) ? d x s u d + s u d-s, 

and 

[ - ^ (t, x) G d+S\d x S 

[%x-vta>> ii(t,x)ed 1 S 

The definition of Qs|?? o (C) is completely analogous to that of Qs\c°{v) an d uses a 
Markov chain f) t ,to ^ t ^ t\ run in the time reversed direction, whose transition 
probabilities are analogous to (12. 61) . (12. 9p . (12. 101) . Then if Ha (77°) is the product 
geometric distribution on configurations (°, we set 

and obtain the equality for the probability measures on Eg 

Ps = Qs = Qs = Ps- (2.14) 

2.2.4 Consistency 

Let S' C S" be two bounded hexagonal domains in Z 2 , and let Ps>,Ps>> be the 
probability distributions on the configuration spaces Eg', Eg//, respectively, as de- 
fined above. The probability measure Ps>> on Eg" induces a probability measure 
Ps"\s' on Eg' which is the -Pg» — distribution of the restricted process fj tjX , (t, x) G S. 
Then the following consistency property is true: 

Ps»\s> = Ps>- (2.15) 

The proof of (12.151) uses (12.141) and the argument in Arak and Surgailis |AS89al 
Theorem 4.1]. Let Xh X/i X\ denote a vertical line, an ascending line, a descend- 
ing line, respectively, in Z 2 , with the last two having slopes 7r/4, — 7r/4, respec- 
tively. Any such line x partitions Z 2 into the left part Z 2 _ = {(£, x) G Z 2 : t ^ 
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t', x ^ x' for some (tf, x') G x} (which contains the line x itself) and the right 
part Z 2 )+ = J?\L\_. Note it suffices to show (l2~T5l for 

s' = s"nz 2 Xt± , 

for any line x = X/> X\ of the above type. 

In the simplest case S' = S" H _ (12.151) trivially follows from the construction: 
in this case, Ps"\s' is nothing else but the evolution 77*, to ^ t ^ t' observed 
up to the moment t' ^ ti,(x,t') G Xh anc ^ therefore coincides with Pgi. The 
case S" = S" fl 2^ follows by the observation that 77^, 6 5(1 Z^ />+ do 
not participate in the definition of the probability of fjt, x , (t, x) G S fl Z 2 _: the 
evolution of the particles after they exit through x / h &s 110 effect on the evolution 
before they exit this line. The case S' = S" fl Z£ _ is completely analogous. 

Non-trivial cases of ( 12. 15ft are 5" = S"' fl Z^ + , x — Xt^X/jXX- However, the 
'reversibility' (12. 14ft allows to exchange the right and left directions by replacing 
Ps by the 'reversed' process Ps and thus reducing the problem to the above 
considered cases of S' = S" fl Z£ ■ 

By consistency f !2. 15j) . the evolution of particles defined in finite hexagonal do- 
mains, can be extended to the evolution on the whole lattice Z 2 . Let £ = 
be the set of all configurations r\ = {j]f x) r]^ x , (t, x) G Z 2 } satisfying (12. 4ft for 
each (t,x) G Z 2 . Then there exists a probability measure P = Pfp. on S whose 
restriction P\s to an arbitrary hexagonal domain S coincides with Pg\ 

P\s = Ps- 

Furthermore, P is invariant with respect to translations of Z 2 . 
2.2.5 The discrete broken line process 

In this section we shall describe the construction of the broken line process in 
a finite hexagonal domain S C Z 2 . This construction is more similar to the 
construction in the continuous Poisson model, and will yield a more streamlined 
proof of the the law of large numbers of the intersection process defined below. 

Given a hexagonal domain S C Z 2 , we denote by S = Ut, x) G Z 2 : \(t,x) — 
(t',x')\ < y/2 for some (t',x') G S 1 }. In other words, S consists of all points of S 
plus the neighboring points which lie at the distance \/2 to the boundary of S. A 
broken line in S is a finite path 7 = {(tj, i) : (tj, i) G S} such that — tj| = 1 
for all i G Z such that (ti,i) G S, (ti + \,i + 1) G S. A broken line in S can be 
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identified with the graph of a continuous function x ^ : ^-(7), h(t)] ~^ K 
which is linear on each interval [z, z + 1] and satisfies 7(2) = £j for each i e Z; here 
i_(7) < ^+(7) are integers which satisfy (tj ± ( 7 ), ^±(7)) G ^S 1 . Clearly, the graph 
of the function 7(-) belongs to the union of all edges of S which are identified 
with closed intervals of the length \f2 between nearest neighbors of S. For each 
(t, x) G S, we shall denote by ef tX the edge (interval) between the points (t, x) and 
(t + l,x±l). 

Given two broken lines 7', 7", we write 7' ^ 7" if ^'(x) ^ 7"(x)Vx G [z_(7'), 
z+(V)] D [i_(7"),H(7")], and 7' -< 7" if 7' ^ 7", 7' ^ 7" holds. For any (t, x) G 5, 
we denote by ef x ,e^.,e^ x ,ef x the edge of Z 2 incident with (t,x) and lying in 
northeast, northwest, southeast, southwest direction, respectively, so that e{ x = 
x+ i, e^. = e^j x+1 . Clearly, each edge incident to a point of the form (t, x) G 
S, connects this point with some point of S. 

With each configuration fj = ((°,r]) G (see Section l2.2.2p . one can associate a 
finite partially ordered system {7^, j = 1, . . . , M} of broken lines in S, such that 
for any point (t,x) G Z 2 , the relations 

v(eQ = vtxi K e i) = VT,x, ^( e i) = C^, ^(eC) = Ci> ( 2 - 16 ) 

hold, where z/(e) is the number of broken lines 7« which pass through a given 
edge e of Z 2 , and where r]f x1 r]j~ x , (f x , (f x are related by fl2.ip . fl2.3p . and denote the 
respective numbers of outgoing and incoming particles to a given site (t, x) G S, 
see Section 12.2.11 Because of these relations and (12.161) . it is convenient to denote 
configurations fj G Ss as 

V={ V (e),eeS(S)), (2.17) 

where S(S) is the set of all edges of S, and where 77(e) is the number of 
particles which pass through along the edge e; 77(e) = vtxiVtxiCtxiCtx f° r 
e = ef x , ef x , ef x , e^ x , respectively. This field 77 is what we shall call flow field. 

According to this notation (12.161) becomes 

v(e)=7i{e), ee£(S). 

The broken lines for a given configuration fj G £s can be constructed as follows. 
For each particle passing through a given edge e G E(S), we define a label p G 
{1, . . . , 77(e)}, which is interpreted as the relative age of that particle among all 
particles which pass through the same edge. The particle whose relative age is 
p — 1 is the oldest and that with p = 77(e) is the youngest. This way, we put a 
coherent order on moving particles; more formally, any particle is characterized 
by a pair 

(e,p), e G E(S), p= 1,2,..., 77(e). 
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We shall define now a relation between labeled particles (e\,pi), (e2,P2) on two 
adjacent edges e%, e 2 G E(S), which we denote by 

(ei,Pi) ~ (e 2 ,pz). 

If relation (ei,pi) ~ ( 62,^2) holds, we say that they are associated, or belong to 
the same generation. The adjacent edges ei,e 2 by definition are any two edges 
e t,x> e t^r> e <vr> e C incident with some (t,x) G S. The relation (ei,pi) ~ (^2,^2) is 
defined in the following cases: 

Case 1: e x = e( x , e 2 = e^ x ; p 1} p 2 < r/( e i) A r?(e 2 ) 

Case 2: e x = e^ x , e 2 = e t x x ; r/(e 2 ) > p 2 > r/(e(l). 

Case 3: e l = e( x , e 2 = e{ x ; 7]{e x ) > rj(e^ x ), p 1 > r](e^ x ). 

Case 4: e x = e^ x} e 2 = e( x ; rj(ei) - £ t ,x < Pi < v( e i)- 

Namely, (ei,pi) ~ (e 2 ,p 2 ) holds if and only if 

Case 1: p 2 — Pi] 

Case 2: p x = p 2 - r](e( x ); 

Case 3: p 2 = pi -r){e^ x )\ 

Case 4: ^(ei) — p% = r](e 2 ) - p 2 - 

It is clear that associated particles are either the particles which annihilate each 
other at (t,x) (Case 1), or the particles born at (t,x) and moving in opposite 
directions (Case 4), or, as in Cases 2 and 3, we associate a younger incoming 
particle which is not killed at (t, x) with the corresponding older outgoing particle 
(both moving in the same direction). 

Now, we say that a finite path 7 = {(tj, i) : (t i; i) G S} is a broken line in a given 
configuration fj (12.171) if any two adjacent edges of this path are associated, i.e., 
belong to the same generation. Let {jj, j = 1, . . . , M} be the family of all broken 
lines in a given configuration fj. The family = 1, . . . , M} is partially order 

by the relation ^ defined above (with the possible exception of broken lines that 
exit S at its left or right boundaries, in which case two lines may not be ordered). 
This follows from the definition of a broken line and the fact that different pairs 
of associated particles cannot 'cross' each other: if (e'^p'J ~ (e' 2 ,P 2 ), ( e i>Pi) ~ 
{e 2 ,p 2 ) and e[, e 2 , e'{, e 2 are all incident with the same vertex (t, x), then the paths 
{e'i,e' 2 ) and {ei,e 2 ) cannot cross each other, in other words, they cannot lie on 
different lines intersecting at (t,x). 
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2.2.6 The intersection process 

Let x be an arbitrary line in Z 2 of the type x — X/i X\j Xh X->j where x^ 
stands for a horizontal line: x^ = G Z 2 : x = x'} for some x' G Z. 

For any such line x an d a finite hexagonal domain S* C Z 2 we shall define the 
intersection process as the number of broken lines ji's which intersect x & t some 
point (t, x) G x H 5 1 . However because of the peculiarity of the discrete situation, 
the definition of the intersection process requires some care. We shall consider 
separately the four above types of x- 

1 . Case x — X/ ■ Let 

nAx) ■= K e ^)> (t, x) G x/ n S 
be the number of broken lines which intersect (t, x) 6 X/ from northwest. 

2. Case x = XV Let 

nAx) ■= v(e( x )> (*> x ) G x\ n S 

be the number of broken lines which intersect (t,x) G X\ from southwest; we 
recall that e^ x = e{_ X x _ x is the edge between (t, x) and (t — 1, x — 1). 

3. Case X = XT- Let 

^(x) := K<£x)> r+.(x) := (t, x) G Xt n 5 

be the number of broken lines which intersect X] a ^ {t, x) from northwest and 
southwest, respectively. 

4. Case x = X->- Let 

n + x(x) := K<1), r M (x) := ^(e^), (t,x) G x^ H 5 

be the number of broken lines which intersect x~> a t (t, x) from northeast and 
northwest, respectively. 

The intersection processes are correspondingly defined as 

r S {x/) = (Mx/)) (M)ex/nS 
t 5 (xt) = (^i(XT)^ M (XT)) (M)exT ns 
r s(X\) = (n,x(x\)) ( t ,x)e X \,ns 
t S (x^) = { T tx(X->),T^ x {x->)) {tjX)Ex ^ nS 

The properties below follow from definition (12.161) and the consistency (12. 15ft . 
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The intersection processes Ts(xy), T s(x\)i T s(xi) are i-i-d- sequences of Geom(A)- 
distributed random variables. The process Ts{xS) is stationary, in the sense that 
there exists a stationary process 

T (x^) = { T t^ x (x^)i T t, x (x^)) (t,x)ex-> J 

such that, for each finite hexagonal domain S C Z 2 , the distribution of the process 
t(x~>) restricted to (t,x) G S coincides with the distribution of Ts{xS)- 

The marginal probability laws of all four processes coincide, i.e. 

Tt,x{X/) = n,x(x\) = ^(Xt) = r i(^-) = Geom(A). 

Note that the process t(x~^ is not i.i.d. This fact follows by an explicit compu- 
tation of the bivariate distribution 

PlT^lxJ) = m, r+ (x^) = n) = (1 - A)(l - A 2 )A 2 ( m+ ")(A- m + A"" - 1 - A). 
However, 

£r+(^)=£r t -(^) = A/(l-A). 

Let 

n 

Tn:=J2( T Ux^+r t - (x^)) 
t=i 

be the number of broken lines which intersect interval (0, n) of the time axis (the 
sum above is taken over all points (t, 0) G Z 2 , < t < n). Then a.s. 

lim n~ l T n = A/(l - A), 

n— +oo 

that is, the intersection process satisfies a law of large numbers. 

2.3 The continuum broken line process on a dis- 
crete domain 

In this section we present a fairly natural generalization of the broken line process, 
namely the continuum broken line process. Within this framework, instead of 
having a certain number £ G Z + of pairs of particles being born at each site, 
we have a mass £ G R+. In this case a broken line is described not only by the 
sites it occupies, this object has some thickness as well. The continuum broken 
lines are not determined by sequentially associating pairs (e,p), we rather need 
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to associate mass. Accordingly, neither can they be understood by just following 
the mass associations over and over, as this mass may well branch due to the 
association rules. At first we consider only broken lines of finite length. But this 
restriction is in fact empty: all the important features of the process are well 
captured by describing finite broken lines. 

Applications to last passage percolation are shown in Section 12.41 In view of 
Proposition 12.11 and its explicit construction, it will hopefully be clear that the 
broken line process and its continuum version are indeed natural objects. The 
proofs of Theorem 12.21 and Theorem 12.31 illustrate how this process can be useful. 

We have chosen to frequently abuse the notation to avoid obscuring the argument. 
Such abuse should not give rise to any confusion. 

2.3.1 Evolution of a continuum-particle system on the dis- 
crete lattice, flow fields 

Let S be a fixed hexagonal domain in Z 2 . Consider S defined on p. ESI and take 
dS = S\S, £(S) = {e = (y,y') e £{I?) :yeS). Usually we will denote the 
elements of Z 2 by y = (t, x) and the elements of £(Z 2 ) by e. We call S a rectangular 
domain if it is a degenerate hexagonal domain, i.e., #<9o£ = #d%S = 1. It is 
convenient to define, for a rectangular domain S, the boundaries d±S = {(t,x) G 
dS :(t + l,x^l) eS} and d ± S = {(t,x) 6 dS : (t - 1, x =f 1) e S}. 

Consider {C y } y es, £y ^ Vy G S, the particle birth process in S. (° is the boundary 
condition, or the particle flow entering S, as defined on p. [321 One takes r\ y ^ 
0, y e S, the particle flow inside S and exiting S, defined by (12.21) . ( J2~3l) . 

Define fj = {77(e) ) 0, e 6 £(S)}, called the flow field in £(S) associated with 
(C°)0- A s i n the discrete case, there is a 1-1 correspondence between (C°,0> 
((°,r)), and fj. We shall write /7(C°>£) to denote the flow field f] corresponding 
to the birth process £ and the flow £° entering the boundary of S, analogous for 

In general, any nonnegative field r\ that satisfies 

n + — n~ = m + — m~ (2-18) 

for each vertex y, where n + = r](e^), n~ = 77 (e^), m + = r)(e^") and m~ = r)(e^), 
is a flow field. Such a flow field may be defined either on a hexagonal domain 
£(S) or on all £(Z 2 ). In the former case there always exists a unique pair (C°,£) 
such that fj = fj((°,^). In the later case such an expression does not make sense, 
although it is possible to determine £(77) for a given flow field fj. In Section T2.3.3I 
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we shall discuss the extension of rj to all of £(Z 2 ) in the same spirit as for the 
geometric broken line process. 

2.3.2 The continuum broken line process 

Let S be a fixed rectangular domain in Z 2 . For a given flow field fj in S(S), 
we shall define a process of lines whose elementary constituents, called atoms, 
are given by a point p standing at a edge e, i.e., they are pairs of the form 
(e,p),ee£(S),pe (0,77(e)]. 

We associate two atoms and write (ei,pi) ~ (e2,P2) according to the following 
rules: 

Case 1: e x = e{ , e 2 = e^; p x G (0, ^(ei) A r](e 2 )] , p 2 = Pi] 

Case 2: d = e^, e 2 = e^; p 2 G (r}(e{),r}(e 2 )], p x =p 2 -r](e{); 

Case 3: e l = e{ , e 2 = e{] p x G (^(e^), r){ei)\ , p 2 = Pi - v(e^); 

Case 4: e x = e^, e 2 = e^; p< G (r?(e;) - £„, ^(e*)] , r?(ei) - pi = r?(e 2 ) - p 2 ; 

Notice that, given a vertex y, each atom standing at an edge incident to y from 
above is associated with exactly one atom standing at another edge incident to y 
from below and vice versa. Notice also that '~' is not transitive. 

By J we will always denote an interval of the form (p,p'\ and \J\ = p' — p. When 
p ^ p', J = and | J] = 0. We associate two intervals of atoms standing at 
adjacent edges and write (ei, J\) ~ (e 2 , J 2 ) if for all p\ G J\ there is p 2 G J 2 such 
that (ei,pi) ~ (e 2 ,p 2 ) and vice versa. Notice that in this case | J±\ = \ J 2 \. 

We define a broken line as a tuple 7 = (2/0, e±, Ji, yi, e 2 , J 2 , y 2 , . . . , e n , J n , y n ), n ^ 
1, such that 

e» = (Vi-iiVi) 

x i = x i - 1 + l , (2.19) 

ti — ti-1 ± 1 

and yi = (ti,Xi) for i — 1, . . . , n, and such that | J±\ = \ J 2 \ = ■ ■ ■ = \ J n \. If Jj = 
we identify 7 = 0. 

We define next the weight of a broken line w(j). If 7 7^ 0, we put w{j) = \ Ji\ > 0, 
otherwise we let 1^(7) = 0. 

The domain of 7, D(j), is given by ^(7) = {:r , . . . , £„}. Since for each x G .0(7) 
there is a unique t such that (t, a;) G 7 we can denote such t by t(x). (Notice that 
we abuse the symbol 'G' here, since 7 is a tuple and not a set.) 

For 7 = (?/o, ei, Ji, yi, e 2 , J 2 , y 2 , . . . , e n , J n , y n ) we define what is called its trace by 
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%) = (l/o, ei, 2/1, e 2 , 2/2, • • • , e n , y n ). In general, any tuple £ satisfying (|2.19|) will 
be called a broken trace. Since either (y , . . . , j/ n ) or (ei, . . . , e n ) are sufficient to 
determine £, we shall refer to any of these representations without distinction. 

The domain of a broken trace £ is defined as D(£) = {x : (t, x) G £}. (Again 
we abuse the symbol '€'.) We also define I(£) = {t(x) : x G D(£)}. It follows 
from (12.191) that D(£) is convex, i.e., it is of the form {a, a + 1, a + 1, . . . , b — 1, b} 
and that the same holds for I{£). 

We write £ C S if y e SUd^SUd'S, y n e SUd+SU d + S and y u . . . , y n _i G 5, 
or, equivalently, if ei,...,e n G £(5). We say that £ crosses S if £ C. S and 
2/o, 2/n e 95. Let C(S) = {£ C S : £ crosses 5}. We write 7 C 5 if %) C 5 and 
we say that 7 crosses S if £(7) crosses 5. 

We say that the broken line 7 C 5 is associated with a flow field fj defined on £(S) 
if (e,_i, Jj-i) ~ (e*, Jj) for i = 1, . . . , n. B(fj) will denote the set of all broken lines 
associated with fj. 

For V C Z 2 , and £ = (2/0, . . . , y n ), define £ HV = {y { : y, e V,i = 0, . . . ,n}. 

For a given broken line 7 associated with the field 7/(C°,0> its left corners (£, x) 
correspond to part of the particle birth £ tjX at (t,x). Given a broken trace £, we 
denote by L(£) the set of left corners of £, i.e., the points (t,x) G £ such that 
(t + l,x±l) e £. Also let L( 7 ) = L(%)). 

We define the fields and £(7) in Z 2 by 

m]v = *W)(v) (2-20) 

and £(7) = w(7)n L(7) . 

In the same fashion, the extremal points of a broken line that crosses S correspond 
to a particle flow entering or exiting S. So we also define 



1, (x,t) G £f] S and (t - l,x- 1) G £ H <9S 
otherwise 



[C" 



(x, t) G £ n S and (t - 1, x + 1) G I D 95 
otherwise 



take the corresponding £° and denote by (°(£). Define C°(l) = w (l)(°(^(l))- 
Also define fj(£) by fj(£)(e) = H ee ^ and 77(7) = u> (7)^(^(7)). It is easy to see that 
fj(£) = f/(( (£)^(£)) for £ G C(5), analogous for 7. 

Let 77 be given and fix some £ = (y , e±, 2/1, e 2 , 2/2, • • • , e n , 2/n) ^ 5. There is one 
maximal broken line that has trace £ and is associated with the field fj, which will 
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be denoted 7(^), the dependence on 7] is omitted. By this we mean that there 
exist unique J 1; J 2 , . . . , J n such that j(£) = (y , e u J u y u . . . , e n , J n , y n ) G B(ff) 
and such that any pi,p 2 , • • • ,p n with property (ei,pi) ~ (e 2 ,p 2 ) ~ • • • ~ (e n ,p„) 
must satisfy p, G Jj, i = 1, . . . , n. The proof of this fact is shown in Appendix I A. 1[ 

Let w(£) denote the maximum weight of a broken line in B{ff) that has trace £, 
which is given by wfj (£)'). The dependence on the field fj is omitted when it is 
clear which field is being considered, otherwise we shall write w v (£). Notice that 
D ("((£)) = D(£) if w n {£) > 0, and D(j(£)) = otherwise. 

We write £ C £' if £' — (y , ei, yi, e 2 , 2/2, ■ ■ ■ , e n , 2/n) and, for some ^ a < b ^ n, 
£ = (ya,e a+1 ,y a+1 ,e a+2 ,y a +2, ■ ■ ■ ,e b ,y b ); or, equivalently, if D(£) C L>(f) and 
£(x) = £'(z) for all xeD(£). 

Notice that w(t\) ^ w(£ 2 ) when ^ C £ 2 . This is due to the successive branching 
caused by birth/collision that makes longer lines become thinner. It is even pos- 
sible that for £\ C £ 2 C • • • C £ n C • • • we have u>(€ n ) [ and the limiting object 
could have one atom but null weight. 

For broken traces £,£', we write £ y £' if £ is to the right of £'. This means that 
t(x) ^ t'(x) for all x G D(£) n D(f ) and that ^ for some x e a/ G 
D{£'). (The last condition makes sense in case -D(^) H -D(T) = 0.) In general this 
relation is neither antisymmetric nor transitive. Write £ y £' if £ y £' and £ 7^ £'. 

Lemma 2.1. If S is a rectangular domain, the following assertions hold: 

1. Relation y, restricted to {£ G C(S)} , is a partial order. 

2. The elements of C(S) are extremal in the following sense. If £ G C(S), 
£' CS and£C£>, then £ = £' . 

3. Relations C / y have a certain concavity in the following sense. Suppose 
that £ C h and £ C £ 3 /or some £ C S and £ x < £ 2 r< 4 e C(S); toen £ C £ 2 . 

It should be clear that the items above hold true. The rigorous proof of this 
lemma consists of straightforward but tedious verifications and is postponed until 
Appendix IA.21 

Lemma 2.2. Let S be a rectangular domain. If w n (£) > and w n (t') > for 
some flow field fj, then £ and £' are comparable, that is, £ y £' or £' y £. 

Lemma [2T21 is a consequence of the way the association rules have been defined. In 
order to prove it, one keeps applying such association rules and the result follows 
by induction - see Appendix IA.31 
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The theorem below is fundamental. It says we can decompose a given flow field in 
many other smaller fields, each one corresponding to one of the maximal broken 
lines that cross S and is associated to the field. In this case the original flow field 
and all of its features, namely the birth process, the boundary conditions and the 
weight it attributes to broken traces, are additive in the sense that each of these 
is obtained by summing over the smaller fields. On the other hand it tells that, 
given an ordered set of broken lines, it is possible to combine the corresponding 
fields and thereby determine the flow field that is associated to them. 

Theorem 2.1. Let S be a rectangular domain. 

Given l\ -< ■ ■ ■ -< £ n £ C(S) and wi, . . . , w n > 0, there is a unique flow field fj in 
S(S) such that 

= %S i = l > S ° TSOmei foranyeeC(S). (2.21) 

I (J, otherwise 

Moreover, for such that fj = fj((°,£,), the following decompositions hold: 

'/ E "yir C° = 5>iC°, £ E "'A,- (2-22) 

j 3 j 

where fjj = fj(£j), Q = C°(^j) an d Cj = Furthermore, 

w v (£)= Wr,{£')-a.ici' foranyiCS. (2.23) 

t'ec{s) 



Conversely, given a flow field fj, there are unique sets {£j} and {wj > 0} that 
satisfy 112. 21\) . The set {£j} is totally ordered by the relation V ' and \2.22\) - [2~2^\) 
hold in this case. 



Proof. We start by the converse part, first proving that (I2.23f) holds for any 
flow field fj, which is the most laborious work. The proof of (I2.23P is a mere 
formalization of the construction described below, whereas (12.211) and ( 12. 22ft are 
immediate consequences, as discussed afterwards. 

Then it will suffice to show that, given any pair of sets {£± -< i% -< ■ ■ • -< £m} 
and {wi,...,wm > 0}, there is some flow field fj satisfying ( 12.211) . Unique- 
ness of such flow field follows from the converse part. Again, by the converse 
part, ( 12.221) and ( 12.231) will hold as well, completing the proof of the theorem. In 
order to prove that (12.211) holds we basically have to see that the same construction 
can be reversed. 
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Though the construction looks simple, several equivalent representations of a flow 
field may be seen on the same picture and the theorem will be deduced from these 
representations. We remind that the sole condition for a field 77(e) ^ 0, e G £{S) 
to be a flow field is the conservation law below: 

V&) + V(e{) = V(e{ ) + V(e}) V y G S. (2.24) 

The construction we start describing now strongly relies on this fact. 

First we plot on a horizontal line two adjacent intervals whose lengths correspond 
to the flow 77 on the two topmost edges of S(S), i.e., the two edges incident from 
above to the topmost site y G S. Since the intervals are adjacent we may do that 
by marking three points on this line. On the next line (parallel to and below the 
first one), we plot four adjacent intervals having lengths corresponding to the flow 
on the four edges incident from above to the two sites on the 2nd row of S. It 
follows from (I2.24p that one can position these intervals in a way that its 2nd and 
4th points stay exactly below the 1st and 3rd points of the first line. By linking 
these two pairs of points we get a "brick" that corresponds to the topmost site 
y of S. The width of this brick is equal to the quantity expressed in f)2.24p . its 
top face is divided into two subintervals having lengths 77 (e^) and fj{e^) and its 
bottom face is divided into intervals of lengths r](e^) and //(e^) - see Figure I2~2l 

We carry on this procedure for the 3rd horizontal line, thereby getting two bricks 
that correspond to the two sites on the second row of S. We keep doing this 
construction until we mark intervals corresponding to the two bottommost edges 
of S(S). In the final picture we have one brick corresponding to each site of 5*. 
Each pair of bricks on consecutive levels that have a common (perhaps degenerate) 
interval on their boundaries corresponds to adjacent sites y and y' in 5"; the length 
of this common interval equals 77 ((y, y')). Once all the intervals have been plotted 
at the appropriate position, forming all the bricks, we draw a dotted vertical line 
passing by each point that was delimiting these intervals. By doing this we divide 
the whole diagram into strips, completing the construction of the brick diagram 
as shown in Figure 12. 2[ 

Each strip corresponds to a maximal broken line that crosses S and the weight of 
this broken line equals the width of the strip. The sites/edges that compound the 
broken line correspond to the bricks/intervals the strip intersects. Given a broken 
trace £ C S, we can determine the maximal broken line ^(l) that passes through 
£ by looking which strips pass by all the sites/edges of £ (i.e., their corresponding 
bricks/intervals); the weight of this broken line is obtained by summing the width 
of such strips. See Figure 12.31 

So, given a flow field fj we construct the brick diagram from which the broken 
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Figure 2.2: construction of the brick diagram. In this example S is a 3 x 3 rectangular 
domain. There are 6 horizontal lines with a total of 24 intervals forming 9 bricks. The 
diagram is divided into 15 strips. 

line configuration is deduced, with the desired property that such broken lines 
satisfy (12T2TD and (12T23D . 

On the other hand, let a set of well ordered broken lines, that is, {i\ -< £2 ~< ■ ■ ■ -< 
£ M } and {wi, . . . , % > 0}, be given. One can consider the corresponding broken 
line diagram, from which one constructs the brick diagram and the later gives a 
flow field satisfying (12.211) . 

In a first reading, one is encouraged to understand the above description with the 
corresponding figures. The more interested reader will find a detailed proof in 
Appendix IA.41 

Since (12.231) has been proven for any flow filed, notice that uniqueness of {£j} and 
{wj} is trivial by definition and well ordering follows from Lemma 12.21 Finally, 
( I2.22j) holds because of ( 12.23!) when we write each process in terms of weights 
of certain broken lines. For C°,^,r] such that fj = fj((°,l;) = fj((°,i]) we have 
Cy = V(e^), = v(e{), % = V(e^), Vy = v(e{), and £ y = T]+ A 77" . Given 
e = {y,y') E S(S), take £(e) = (y,e,y'). Then of course 77(e) = w v (£(e)). Also, 
given y e S, take £<(y) = e_, y, e+, y+), where y = (t,x), y± = (t + 1, x ± 1) 
and e± = (y,y±). It is also easy to see that w ri ( K £ < {y)) = £ y . Also notice that 
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Figure 2.3: broken lines that cross S and calculation of the weight of a given broken 
trace. On the left we have a random configuration of broken lines on the 3x3 domain 
S and on the right we have the brick diagram from which they arose. The 3rd strip 
in the brick diagram is highlighted and the corresponding broken line appears in bold 
on the left. Also, for a broken trace £ that starts at the center of S, goes northeast, 
and then exits S going northeast again, we have determined which broken lines in C(S) 
contain £; they are the 11th and 12th ones and they also appear in bold. In this case 
w(£) equals the sum of the widths of the 11th and 12th strips. 



( y (£) = Hg"^ = ll^ e ^ cr Now we put all the pieces together to get ( y = 

rj(e^) = w„(*(e\)) = Y,jU } j\$>)ct j = Ej w AQ ).</• Le -> C" = Ej w d Q ■ Th e 
other equalities are deduced similarly. 

This completes the proof of the converse part and, as discussed above, of the 
theorem. □ 



Corollary 2.1. Let S be a rectangular domain and let r]((°,£) be given. Then 



E C+ E 

yed-S yed'S 



E c + E < 

y£d+S yed+S 



E 



(2.25) 



Proof: For each £ G C(S) there is exactly one e G £ such that y G d + S, y' G d + S 
or y G d + S,y' G 9+5, where e = (y,y'). So 1 = E^es+s R((t,x),(t-x,x+i))ee + 
E(t,*)e0+s %t,*),(t+i.*+i)>e< = E„ ea+ sCW + E^e+s^W- Multiplying by w{£) 
and summing over all £ G C(S) we get J2 y ed+sCy + J2 y& a+ s Vy = E te c(s) w W- 
The proof of Eyea-S + J2 ye a-sVy = E te c(s) W W is similar. □ 
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Denote by Hs((°,t;) the quantity expressed in ( 12.25ft . if C° = we just write 

For a birth field £ ^ on S we define the last passage percolation value G\s(£) as 
the maximum of ^2 yen £ y over all 7r G lis, where 

lis = {n = (yo, ■ ■ ■ , y m ) ■ yo e doS 1 , y m e dxS, t i+1 = ti + i, x i+1 = ^ ± l}. 

When it is clear which rectangular domain we are referring to we shall drop the 
subscript S of H s , G$ and rig. 

The next proposition illustrates the connection between broken lines and pas- 
sage time. Furthermore, its proof gives an explicit algorithm for determining the 
optimal path (which is a.s. unique when £ has a continuous distribution). 

Proposition 2.1. Let S be a rectangular domain. The last passage percolation 
value Gg(£) is given by the sum of the weight of the broken lines associated to the 
corresponding birth field: Gs(£) = Hs(£). 

Proof: The proof consists on formalizing the following argument. On one hand, 
an oriented path tt connecting 8qS — > d±S can cross at most one left corner of each 
broken line. On the other hand, it is possible to assemble the path backwards, 
following a local rule that does not miss any broken line. This is possible because 
the broken lines never cross each other. See Figure l2~4l 




Figure 2.4: construction of the maximal path. From £ one constructs the flow field 
and by following algorithm (|2.27p one gets an optimal path. The theory developed for 
broken lines and Theorem 12.11 guarantee this is indeed optimal. The algorithm forbids 
the path to cross any broken line; for that purpose it suffices to require that it does not 
'cross' the flow field. 



Take fj = fj{(° = 0,0 and write {£ G C(S) : w v (£) > 0} = {h -< £ 2 -< < M- 
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Now by ( 12.201) and (12.221) it follows that G(£) is the maximum over all n G II of 

Since the paths ir are oriented, they cannot intersect more than one left corner 
of each broken line, hence ^2 yG7r ^-L(ej){y) ^ 1 for each j. We shall exhibit an 
algorithm for constructing a path 7r* that satisfies 

X>J*W)(V)>1> (2-26) 

■(/eyr* 

which completes the proof. 

The path 7r* is constructed by the following rule. Let y m G diS. For % = m,m — 
1, m — 2 . . . , 2, 1, let = t j — 1 and 

x i _ 1 =/ X '- 1 ' rt^^)' (2.27) 
I Xj + 1, otherwise. 

It remains to show (I2.26p . i.e., that tt* intersects L(£j) for each j — 1, . . . , M. Fix 
j and write £ for £j. Assume without proof that tt* intersects £ at some point (a 
complete proof is shown in Appendix IA.5H . Take n = min{i : G £}■ If yo G ^ we 
have E L(£) due to the fact that (y = 0. So suppose n > 0. By construction 
x n-i — s n ±l; assume for simplicity = x n + 1. Now t n _i = t n — 1 and 
n is minimal, so (t n — l,x n + 1) £ Since y n g" dS, x n + 1 G -D(£) and thus 
(t n + l,x n + l) G £. Now if (t n — l,x n — 1) were in £ there would be p\ G (0, T]{ e y n )]i 
P2 G (0, 7/(e^)] associated by Case 3, which implies v( e ^ n ) > r /( e ^)' contradicting 
the choice of x n _i = x n + 1, so (t n — 1, x n — 1) ^ £. But also x n — 1 G -D(£), thus 
(t n + 1, x„ — 1) G £ and therefore y n G L(£). □ 



Corollary 2.2. H((°,^) in nondecreasing in (C°>£)- 

There is a nice direct proof for the above corollary that uses Theorem 12.11 and a 
much shorter one that relies on Proposition EHJ The argument for the short proof 
will appear in the proof of Theorem 12.21 just before (12.351) . so it is not repeat it 
here. The next corollary is indeed an immediate consequence of Proposition 12.11 
although not at all obvious when one thinks of the definition of H${£) by itself. 

Corollary 2.3. Let S be a rectangular domain and (^0 o particle birth process. 
Consider the map M : (t,x) i— > (—t,x), S' = M(S), £,'(t,x) = £(—t,x). Then 
H s (0 = H s ,(?). 
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2.3.3 Reversibility 

In analogy with the discrete case, given the distributions of ( + , (~ and £ satisfying 
( 12.131) . the evolution of the particles described as the evolution of the associated 
Markov chain is reversible. More precisely, reversibility is determined by the 
following relation 



In general, all we have to check is that the relation (I2.28P holds for the distribution 
q(dn + , dn~ \m + , m~), where the q is the kernel that makes (12.181) happen a.s. and 
(n + A n~) be independent of m + and vrT . (It is given by £ y , having law 7r 3 .) 

Thus any triple 7i"i,7r2, 7T3 that satisfies (I2.28f) will define a family of measures 
{Ps : S hexagonal domain} which is consistent, i.e., satisfies (12.151) . In this case fj 
can be consistently extended to £(Z 2 ). As a consequence, if one takes {(y) y ed-S 
i.i.d. distributed as 7Ti, (Cy) y ed+s i-i-d. having law 7r 2 , and (C, y )y£s i-i-d. with law 
7r 3 , then {rjy) y £d+s will be i.i.d. with law 7Ti and ij]~) y( zQ-s will be i.i.d. with law 

The geometric broken line process described in Section 12.21 is a particular case of 
r] + ~ Geom(Ai), rf ~ Geom(A 2 ), £ ~ Geom(AiA 2 ), taking Ai = A 2 = A. Below 
we characterize the distributions that satisfy (12.281) . 

Since the particle system always evolves by keeping the relation (12.181) . we can 
parametrize this hyperplane in by T : — > T(r,s,t) = (r,s,t + [r — 
s} + ,t + [r — s}~). The joint distribution of (n + ,n~ ,m + ,m~) is given by the left 
hand side of (I2.28p . It can be obtained by taking the projection T^/z, where 

fl = 7Ti X 7T 2 X 7T 3 . 

Now consider the operator in given by L(x, y, z, w) = (z, w, x, y). Reversibility 
(I2.28P in our case just means that L preserves T*/i. This is equivalent to the fact 
that 



where R = T~ X LT . Writing R explicitly gives R(r, s,t) — (t + [r — s] + ,t + [r — 
s]~,r A s) and R 2 = I, the identity operator. Note that £>i, dB\, B 2 , given by 
Bi — {(r, s,t) : r > s}, B 2 = {(r, s,t) : r < s}, remain invariant under R and they 
form a partition of IRi. Also note that {r A s = 0} is in bijection with {t = 0}. 

Another example that satisfies (12. 29ft and thus (12.281) is when the 7Tj's are respec- 
tively the exponential distributions exp(aj) with a 3 = cti + a 2 . In this case it is 
a simple exercise to check (12.291) . Assume A C Bi, the expression for R becomes 



7r 1 (dm + )7r 2 (dm )q(dn + ,dn \m + ,m ) = 
= 7r 1 (dn + )7r 2 (dri _ )g(dm + , dm~ \n + , n~). 



(2.28) 



R*fi — /i, 



(2.29) 
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much simpler and it is straightforward that fi(R 1 {A)) = [i(A). For AC5 2 the 
situation is analogous and n(dBi) = 0, completing the proof of (j2.29p . 

The question of finding triples of measures on M + that satisfy f 1 2 . 2 9 1) is an inter- 
esting question on its own. In discussion with A. Ramirez we have shown the 
following. If the measures are supported on Z + and assign positive weight to each 
n G Z+, then they must be geometric distributions. If the measures have continu- 
ous positive densities, then they must be exponential distributions. Of course a lot 
remains between being supported on Z + and having continuous positive density, 
but nonetheless we believe that the only examples are these two, modulo adding 
or multiplying by constants. 



2.4 An application: geometric and exponential 
last passage percolation 

It is easy to see that the last passage percolation model satisfies a law of large 
numbers by super-additivity arguments. However it is interesting that for the 
two-dimensional model and for the special case of i.i.d. geometric or exponential 
passage time distributions there is an explicit expression for the limiting constant. 
For the exponential distribution it was found by Rost [Ros81j and for geometric 
case by Jockusch, Propp, and Shor [JPS98J . Large deviations were studied by 
Johansson [JohOOj (lower tails) and by Seppalainen |Sep98| (upper tails). Fluctu- 
ations were studied in [JohOOj . 

With the aid of the broken line theory developed in the previous sections it is 
possible to re-obtain the explicit constants for the law of large numbers. We 
also prove exponential decay for the probability of deviations. This application is 
based on the law of large numbers for the intersection process, on Proposition 12.11 
and on Corollary 12.21 

In the same spirit, O'Connell |O'C00] also devises such constants by simple proba- 
bilistic arguments. We note that our approach is self-contained, except for using of 
Cramer's theorem for large deviations of i.i.d. sums. A proof of Burke's theorem 
is implicitly contained in our considerations of reversibility. 

The construction consists on first choosing the appropriate distributions of the 
boundary conditions that (i) make the broken line process reversible and (ii) pro- 
vide the correct asymptotic behavior, and then dropping the boundary conditions 
afterwards. 

What we present is an alternative proof that could give some geometric insight to 
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the model. Besides, the broken-line approach provides an explicit, linear algorithm 
for determining the maximal path (see the proof of Proposition 12. ip . Our results 
also show that the boundary conditions give no asymptotic contribution to the 
total flow of broken lines that cross a given domain when they have this suitable 
distribution. 

For N,M G N, we define the last passage percolation value on the square 
{1, . . . , N} x {1, . . . , M} as the random number G(N, M) given by the maxi- 
mum sum of £ yj over all oriented paths (yi, . . . ,i/n+m-i) from (1, 1) to (N,M). 
G(N, M) is random because so are the £ y 's. 

Theorem 2.2. Suppose £ y , y G N 2 , are i.i.d and distributed as exp(a), a > 
and let 3 > be fixed. Then a.s. 

lim i GW L/Wj) = 
For each 5 > 0, there exists c = c(5) > such that 

,2 



P 



a 



> 5 > ^ e 



-cN 



(2.30) 



for all N EN. 



Theorem 2.3. Suppose £ y , y G N 2 , are i.i.d and distributed as Geom(A), A G 
(0, 1) and let (3 > be fixed. Then a.s. 



For each 5 > 0, there exists c = c(S) > such that 



1. 



-G{N, l(3N\) 



(1 + v^A) 



for all N EN. 



> 8 y< e 



-cN 



Proof of Theorem 12.21 The central idea of the proof may be hidden among all 
the calculations, basically it is built on the following argument. 

By Proposition 12. 14 Corollary 12.21 and the laws of large numbers for i.i.d. expo- 
nential r.v.'s, 



G(N,(3N) = Hs(0^H s (C,0= £ tf + E % 



yed-S 



y&d-S 



N BN 

— + —; 
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where S is a rectangular domain with N x (3N sites. Here a+ and a_ can be any 
pair of positive numbers that make a = a+ + a_ and therefore the broken line 
process reversible when the are distributed as exp(a±) and the £ are distributed 
as exp(a). As a consequence we have 



lim ^=G(N, (3N) ^ inf 

N a+,a- 



a+ a_ 



(l + v^) ; 



a 



the infimum being attained for a± = a/(l + /3 ±1//2 ). Now we want to have the 
opposite inequality. We argue that Hs(C°,0 cannot be much bigger than Hs(£). 
To compare both, consider, instead of boundary conditions in dS, a slightly 
enlarged domain S' without boundary conditions but where to each extra site we 
associate £' corresponding to the previous so that Hs>(£') = Hs{C°,0- Now 
for Hs'{C) to be considerably bigger than H s(£), it must be the case that the 
^'-optimal path in S' occupies a positive fraction e of the boundary d + S' and then 
takes a ^-optimal path in the remaining (N x (1 — e)/3iV)-rectangle. But in fact it 
happens that an oriented walker, after visiting e(3N sites in d+S, looks ahead and 
realizes it is far too late to perform the ^'-optimal path, as the following equation 
shows: 

a a a a 

Here 

/(r) = l-r/2-vT=7 

is positive and increasing in (0, oo). Therefore the ^'-optimal path in S' cannot 
stay too long at the boundary of S' and thus #s(£) ~ Hs'(£) = Hs(C°,0- 

Now let us move to the proper mathematical proof. 

We shall use the following basic fact. Given a±, ai2, «3, p, S, K > 0, there exist 
positive constants c, C > such that if (X,), (Yj) and (Zj) are sequences of i.i.d. 
r.v.'s distributed respectively as exp(ai), exp(a2) and exp^), then for all JVGff, 



3/, m, n e {0, . . . , \pN] } : £ L *j + E?=i Yj+ 



<: Ce~ cN , (2.32) 



l m_ JL 

Ot\ OC2 Ct3 



>5N -K 



regardless of the joint distribution of (X, Y, Z). 

We first map our problem in Z 2 to the space Z 2 , where the theory of broken lines 
was developed. We do so by considering the rectangular domain S(N, M) given 
by S{N,M) = {(t,x) E Z 2 : ^ t + x < 2(M - 1), ^ t - x ^ 2{N -I)} and the 
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obvious mapping between S(N, M) and {1, . . . , N} x {1, . . . , M}. We write S for 
S(N, M). 

In order to define a reversible broken line process in S(N, M) with creation £ ~ 
exp(a) we can choose a + and a_ such that a + + a_ = a and let C + ~ ex P( a +), 
C" ~ exp(o:_). Take M = M(N) = [(3N\. Choosing a + = a/ {I + f3 1 ' 2 ) and 
a _ = a / {I + I3' l/2 ) gives 



1 , g _ (1 + v^)- 

a i ct_ a 



(2.33) 



By Proposition 12.11 and Corollary [21 



G(^M) = # s (£Ktf s (c o ,o= EC+ E 



and it follows from (12.32112.331) that 



P { ^G{N, M) > + ^ + 5 J> ^ C e- coJV 



(2.34) 



Now let us prove the lower bound to complete the concentration inequality above. 
Consider S'(N, M) = {(t,x) G Z 2 : -2 < t + x < 2(M-l),-2 sC t-x sC 2(iV-l)} 
and for ^ n < /3iV, take S(N,M-n) = {(t,x) 6 Z 2 : 2n ^ t + x < 2(M-1),0 sC 
t-x sC 2(JV-1)}, S'(N,M-n) = {(t,x) G Z 2 : 2n-2 ^ t + x < 2(M-l),-2 $S 
t - x < 2(iV - 1)}, S(iV - n, M) = {(i, x) G Z 2 : < t + x < 2(M - 1), 2n sC 
t-x < 2(iV-l)}, S'(N-n,M) = {(t,x) G Z 2 : -2 < t + x ^ 2(M-l),2n-2 sC 
t-x ^ 2(iV-l)}. 

For (C°, defined on S take £' on S" given by 



Cii,,+i, HM + ie9.5 

otherwise. 



For (C°,0 defined on S take £' on S" given by the analogous formulae. It is easy 
to see that H s (C,0 = H s >(?) and H§(C,€) = ^(O- 
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The two facts below will be important: 

n 

For some n, H s ,(?) = ^ Cj-ij+i + Hg (N>M _ n) (£) or (2.35) 



j=0 



j=0 

For all n, H s (0 ^ #s (J v,M-n)(0 and 

#s(0 ^ #S(jV-n,M)(0- 



(2.36) 



Given any < 5' < 5, by putting (12.351) and (12.361) together we see that for 



G(iV, M) ^ iV 
to hold, we must have either 

or, for some n G {1, . . . , M — n}, 



-5 



a 



(2.37) 



-5' 



(2.38) 



£^Qj+%iV,M-„)(£)^ 

or, for some n G {1, . . . , N — n}, 



;w?) 3 



-8' 



6' 



(2.39) 
(2.40) 

(2.41) 
(2.42) 



The probability of (I2.38P decays exponentially fast and this can be shown exactly 
as was done for ( 12. 34ft . 

We consider now the other possibility, (I2.39|I2.40I) . The case (12. 41112. 42p is treated 
in a completely analogous way and is thus omitted. Let 5 > be fixed, take 
e = 6a/[2P{l + and 5' = § A [^/(e )]. With this choice of parameters 



a 2 



(2.43) 
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holds for e ^ eo and 

^ M , (2 - 44) 

holds for e ^ e . 

Consider the event that (12.401) happens for some ^ n $C |~eo/3iV] = M . Since 

g- 4 exp(a/(l + V / V^)) ^ follows from (|2T43|) and (12T321) that the probability of 
this event decays exponentially fast in N. 

It remains to consider n ^ M and show that in this case it is the probability 
of (12.391) that decays exponentially fast. Now 

P 1 3n G { Mo , . . . , M} : ^ Q + #s (0 ^ iV ^±^L - jV5 j < 
P jsn G {M , . . . , M} : ^ Q + ^(O ^ ^ - N5> 



}< l — J V - I , -, ,_ ^7T\ 2 



j=0 j=n j=0 

where are distributed as exp(a/[l + (/? — eo/?) ±1/ ' 2 ]) so that the broken line 
process on 5 with £ distributed as exp(a) is reversible, and therefore the 77 are 
also distributed as the ( ± . 

By (12.441) the right-hand side of the inequality in the last line is greater than 



and by (I2.32p the last probability above decays exponentially fast. The proof is 
finished. □ 



Proof of Theorem 12.31 The proof is absolutely identical to that of the pre- 
vious theorem, so we just highlight which equations should be replaced by their 
analogous counterparts. 

In the heuristic part take A + = [(A + y/]3\)/(l + G (A, 1) and A_ = A/A+ G 

(A, 1), so that A = A + A_, the process is reversible for ^ ~ exp(A±) and 

A + , 3 A- _ (1 + V/3A) 2 
1-A + P l-A_ 1-A 
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Instead of (I2.3ip consider 



^— + v i-a -^^T^-'-frxW- 



The proof of 



P <( Ic?(iV, l(3N\) > (1 - 1 + * J < e- cJV 

is analogous to the proof of ( 12. 34ft . 

For the opposite inequality, we take e = 5(l-A)/[2(A/3+</A/3)], 5' = § A[^§/(e )] 
so that instead of (I2.43P and (12.441) the following estimates will hold, respectively 
for e ^ eo and e ^ eo: 



px + ^px 6 



e/3 7a^AT + (/5 - e/5) (/3-^)d-A) — + — t^a — < 



1 — A 

The rest of the proof is the same. □ 



Appendix A 
Technical Proofs 



A.l Existence of the maximal broken line 

Here we prove that given a flow field and a broken trace there is a maximal broken 
line associated to that field and having that trace. 

We claim that there exist unique J\ = (ai, 61], J 2 = (a 2 , b 2 ], • • • , Jn — (o n , b n ] such 
that 

(e x , Ji) ~ (e 2 , J 2 ) ~ •- • ~ (e„, J n ), (A.l) 
and such that any pi,p 2 , . . . ,p n with property 

(ei,Pi) ~ (e 2 ,P2) ~ • • • ~ (e n ,p n ) (A.2) 

must satisfy p^ G Jj, i = 1, ...,n. If Ji 7^ we define 7(7) = 
(y Q , ei, Ji,y!,..., e n , J n , y n ) and by (1A.11) we have 7(f) G ^(77); otherwise 7^) = 0. 

To prove the claim start by observing some consequences of the association rules. 

1. If (ei,pi) ~ (e2,p 2 ) then (ei,pi — 5) ~ (e2,P2 — 5) for some 5 > 0. 

2. If (ei,pk) ~ (62,^2) V A; G N with p^ f p x then p\ \ p 2 and (ei,pi) ~ (e2,p 2 ). 

3. If (e!,pi) ~ (e 2 ,p 2 ) and (e 1; pi) ~ (e 2 ,p 2 ) with pi > p x then pi - p' 2 = p x - p 2 
and (ei, (pi,pi]) ~ (e 2 , (p 2 ,P 2 ])- 

Now let A be the set of pi G (0, ^(ei)] for which it is possible to find p 2 ,P3, ■ ■ ■ ,Pn 
such that ( 1A.2I) holds. Suppose that A 7^ and take ai = inf A, 61 = sup A. 
(When A = we take Jj = and the desired properties hold trivially.) Consider 
a sequence (p^) C A with pj 1 f b\. By Property 2 above we have (ei, 61) ~ • • • ~ 
(e n ,b n ) and b\ G A. It follows from Property 1 that a% < b\ and we can take 
another sequence (p\) C A with b\ > p\ \ ai, (ei,pf) ~ ••• ~ (e n ,p£). By 
Property 3 (ei, (pj, 61]) ~ • • ■ ~ (e n , (p£, &„]) and 6* > pf j a, for i = 1, . . . , n; 
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thus (ei, (ai, bi\) ~ • • • ~ (e n , (a n , b n ]) and 61 — ai = • • • = b n — a n . If it were 
the case that a\ > and ai G A, by Property 1 it would hold that a% — 5 G A 
contradicting ai = inf A. Therefore A = (a\, bi]. Suppose pi, . . . ,p n satisfy ( 1A.2I) ; 
by definition pi G A = (ai,&i] and by Property 3 we have bi — Pi = b\ — p 1 G 
[0,6i — ai) = [0, bi — a^, that is, Pi G {a i: bj\. As a consequence we have that 
J'i C (aj, i = 1, . . . , n whenever (ei, J[) ~ • • • ~ (e n , J^), from which uniqueness 
follows. 



A. 2 Proof of Lemma 2.1 



We start by proving Item [TJ Relation y is obviously reflexive. We now show 
that it is antisymmetric. Assume £ y £' y £, so that we are required to show 
that £ = £'. First suppose D(£) n D(£') ^ and take x G D(£) fl £>(f ). Write 
L>(f) = {x_ n , . . . ,x , • • • ,x m }, D(£') = {x_ n r,...,x ,...,x m r} &nd D (£) n D (£') = 
{x-n, . . . , xq, . . . , Xfh} with Xj+i = Xi + 1 and m = m A m', h = n A n'. Since 
£ y £' y £ we have £(xj) = t'(xj) for z = — n, . . . ,rh. All we need to show is that 
n = n' and m = m'. Suppose m = w! . Then ({'(xm),im) G d + S U d+S 1 ; and since 
t{xfh) — t'(xm) we cannot have m > rh, thus m = fh. By the same argument, if 
m = m we conclude m' = rh, therefore m = m = m! . Similarly we show that 
n = n' = n. It remains to consider the case D(£) fl D(£') ^ 0, which is ruled out 
by the following claim. 

Claim A.l. If£,£' G C{3) and D{£) n D{£') = then I{£) n /(f) = 

We will show that I{£) n /(f) ^ implies D(fj n D(£') ^ 0. Let t G /(f) n /(f) 
and take x" G D(£),x' G -D(f ) such that £(x") = t'(x') = t. Assume for simplicity 
x" < x' and let (t, x) denote the topmost site of S. Since (t, x') G S we have 
a; — t ^ x' — t — 2 and x + t x' + t — 2. Now, after passing by (t, x") when 
going upwards, £ must cross either of the lines {(x, t) : x — t = x — t} or {(x, t) : 
x + t = x + t}, because £ G C(S). After crossing either of these lines there will 
be G £ with (x* — 1) — (t* + 1) = x — t or (x* — 1) + (t* — 1) = x + t , 

respectively. Therefore x* ^ x' — But by (12.191) we have 1*\ ^ |x»— x"|, 

so i' - ^ |x* — x"\. Assuming x* ^ x' (the other possibility trivially implies 
the desired result), one has |x' — x*| ^ \x" — x*| thus x* ^ x 1^ and therefore 
G -D(f). Analogously we show that x ' + *" G -D(f ) and the proof is done. 

Finally let us see that y is transitive. For a given point = (t*,x*) G S", define 
A(£*,x*) = {(t,x) G Z 2 : t — x ^ — x*,t + x ^ + x*} = {(t,x) G Z 2 : 
t ^ t* + |x — x*!} and for £ C S 1 , define ^4(f) = U ye ^A(y). It is immediate that 
^4(?/) C iff y G therefore A{£') C A(fj is equivalent to f C A(£). Now 
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let £ y £' y £" . It follows from these observations and from the claim below that 

£ y £". 

Claim A.2. Let £,£' G C(S). Then £' h £ if and only if £' C A(£). 

We start proving the 'only if part of the claim. Let £' y £ G C(S) . Suppose 
D(£) n D{£') = 0. By Claim GO I{£) n I{£' ) = and, since t'{x') > t(x) for 
some x' G D(£'), x G £>(f), we have t'(x') > t(x) for all x' G x G £>(£). 

Consider the case x > x' V x' G D(£'), x G -D(^); the other situation is analogous. 
Writing -D(£) = {xo < • ■ ■ < x n } and to = ^(xo), we must have (to,xo) G d-S 
or (t ,x ) ^ d~S and the later is ruled out since there is x' G D{£') with x' < 
x,t'(x') > t in 5*. Now as (to,x ) G oLi? we have t + x ^ t + x V (t, x) G 
5 and, as x' < xo,t'(x') > t V x' G -D(i") we have t' — x' > to — x for all 
(t',x') G £'. Therefore £' C A(t ,x ) C Suppose on the other hand that 

n D(£') j£ 0, take x G £>(f) n and write D{£) = {x_ n , . . . , x , . . . , x m }, 
D(£') = {x_ n >, . . . ,x , • • • ,x m /} and £>(£) n -D(f) = {x_^, . . . ,x , . . . ,x„} with 
Xj + i = Xj + 1 and m = m A m', n = n An'. For x G -D(^) H -D(^) we have 
t'(x) ^ t(x) and of course (t'(x),x) G A(t(x),x). Take t m = t(x m ), by definition 
(t m ,x m ) G d + S or (t m ,x m ) G <9 + S\ In the later case it must be that m ^ m', for 
if we suppose that m! ^ m, then as t'(x m ) ^ t(x m ), we must have f(x m ) = t(x m ), 
thus (t'(x m ), x m ) G and m! = m. In the former case we have {(t, x) G S : x ^ 
x m } C y4(t m ,x m ). Therefore (t'(x),x) G for x G {xq, . . . , x m i}. Analogous 
arguments show that (t'(x),x) G A(^) for x G {x_ n /, . . . ,Xo}. 

The 'if part is shorter. Suppose £' C A(£). Take some x' G L>(£'), there is 
(t, x) G £ such that (t'(x'), x') G A(t, x) and thus t'(x') ^ t + |x - x'| > t Now for 
x' G £>(f ) n there is some x G D(£) such that (t'(x'),x') G A(t(x),x). We 

want to show that t'(x') ^ t(x') and it follows from the fact that £ C A~(t, x) for 
any (t, x) G £, where A _ (t*, x*) = {(t, x) G Z 2 : t ^ t* + |x — x*|}. 

Item |2] is easy. If £ = (y a , y a+1 , ...,y b )C£' = (y a >,y a > +1 , ...,y' b ) with a' ^ a < b ^ 
6' for £ G C(S), £' C 5, then a = a', 6 = 6' and therefore £ = £'; for if a' < a we 
would have y a G 5, contradicting £ G C(S), same for 6 > b' . Item [3] is not used in 
this work and we omit its proof. 



A. 3 Proof of Lemma 2.2 



Let £,£' C S such that w{i),w{t') > 0. If t(x) = t'(x) Vx G n D(£') the 
result is trivial. So suppose there is xq such that t(xo) 7^ t'(xo) and assume for 



62 



APPENDIX A. TECHNICAL PROOFS 



simplicity that t(x ) < t'(xo). Write D(i) fl D(£') = {x- m , ■ ■ ■ , X-i, xq, x\, Xh}, 
j(£) = {y 

and 

) = {y-m'' e -m'+l> J—m'+M V-m'+li ' ' ' ' e 0' ^0' 2/o> • • • ) e m'' "^ro'j 2/m'}- 

We want to show that t(xi) ^ t'(xi) for z = 0,1, ... ,n; the proof for i = 

0, — 1, . . . , — fh is analogous. 

We claim that, for i = 0, 1, ... ,n — 1, the following facts hold: t(xi) ^ t'{xi), 
t(xi+i) ^ and, in case t(xi) = t'(xi), t(x i+ i) = t'(x i+ i) we have J i+ i -< J' i+1 , 

1. e., p < p' for all p G J%+i, p' G J| +1 . We proceed by induction. For i = 
the result is obvious because of (I2.19P and t(x ) < t'(x ). Suppose the claim 
is true for % = k — 1. We have three possibilities. Case A: t(xk) = t'(xk) and 
t(xk-i) — t'(xk-i). Case B: t(xk) = t'(xk) and t(xk-\) < t'(xk-\)- Case C: 
t(x/i) < t'(xk)- In Case C the claim holds for i = k for the same reason it holds 
for i = 0. In Case A -< J' k and, because of the association rules, we have 
t(xk+i) < t'(xk+i) or t(xk+i) = t'(xk+i) with J fc+1 -< J' k+1 ; either way case the 
claim holds for i = k. In case B, since t(xk-i) < t'(xk~i), by the association 
rules we cannot have t(xk+i) > t'(xk+i) and if t(xk+i) = t'(x k+ i) we must have 
Jk+i -< J'k+i, so the claim holds. 

The proof is complete. 

A. 4 Formalization of the brick diagram 

In this appendix we formalize the construction of the brick diagram described in 
the proof of Theorem 12.11 Then we show that (12.231) holds for any flow field in 
the converse part of the theorem and that (12.211) holds for some flow field in its 
direct part. 

As mentioned in the proof of the theorem, there are several equivalent represen- 
tations of a flow field and this proof draws on each of them. We describe how to 
obtain from each representation the next one, and we mention the properties of 
each representation that guarantee it is possible to come back to the previous set- 
ting. Statements will be made without proof when their verification is a straight 
forward but tedious standard argument. 

To simplify the presentation, assume that S has the form S = {(t,x) G Z 2 : 
\t\ + \x\ < N} for some N G 2Z+. 
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Let S* = {y G Z 2 : 3y' G S : \y - y'\ = 1}. Notice that S* C (Z 2 )* = Z 2 \Z 2 . 

The first alternative representation of a flow field is p t)X , (t, x) G S*, defined below. 

Take W_ N ^i = 0, for i = -N, . . . , take W { = Wj_i + T)(e^ N+i ) and for i = 
1, . . . , iV + 1 take Wi = Wi- X + r){e{_ X N _ i+x ). 

Let p ,jv+i = W . 

For each x = N, N - 1, . . . , 1, 0, -1, . . . , -JV, — iV - 1, set 



for t = N + 1 - x. 

With successive uses of (12.241) it is not hard to see that p t ,x has the following 
properties: 

P-jv-i,o = 0, 

Pt,x ^ Pt+i,a!±i whenever they belong to S* 
and that it is possible to re-obtain fj from the p(t, x) by 



Now consider the set of all p tyX and reorder it by taking A = {q Q , . . . , qu) = {pt, x '■ 
(t,x) E S*} with q < qi < ■ ■ ■ < q M . Then q = W_ N _i = p_at-i,o = and 
Qm = Wn+i = Pn+i,o- 

For (t,x) G S*, take k(t,x) G {0, . . . , M} as the unique subindex k that satisfies 
Pt,x = Qk- Then 



for t 



Pt,x = Pt+l,x+l - V( e t+l,x) 

(N + 1 - x), -(N + 1 - x) + 2, . . . , (N + 1 - x) - 2 and 



r](eS) = p t ,x-i - Pt-i,x 
rj(e£) = p t+ i, x - p t ,x-i 

V( e y) = Pt,x+1 ~ Pt-l,x 

V(e{) =Pt+i,x -Pt,x+i- 



(A.3) 



k(-N- 1,0) = 
fc(iV+ 1,0) = M 

k(t, x) ^ + 1, x ± 1) whenever they belong to S* 
For all k — 0, . . . , M there is y G 5** such that fc(j/) = fc. 



(A.4) 



Of course it is possible to recover p t>x from A and k(t,x): 



Pt,x = Qk(t,x)- 
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For y = (t, x) G S, take k±(y) G {0, . . . , M} as 



'k(t-l,x), y G SUd+SUd'S 

k-(t,x) = { k(t,x- 1), yed+(S) 

k(t,x + l), yed-(S). 

'k(t + l,x), yeSUd+SUd-S 

k-(t,x) = { k(t,x-l), yed + (S) 

k(t,x + l), yed-(S). 



Then 



fc_(-iV,0) = 
k + (N,0) = M 

k+(t,x) < k+(t + 1, x ± 1) whenever they belong to S 
k-(t, x) ^ k-(t + 1, x ± 1) whenever they belong to 5 

x) = + 2, x) whenever they belong to S 
For all = 1, . . . , M there is y G S such that = k 

For all fc = 0, . . . , M - 1 there is y G S such that = fc. 

To obtain fc from fc ± , take, for (£,x) G 5*, 

'0, t = -N-l 

k(t,x) = { M, t = N + 1 . 

+ 1, x) = k + (t — 1, x), otherwise 

We may also consider t(j, x), j G {1, . . . , M }, x G {— iV, . . . , N}. For such (j, x), 
take 

{t , (i, x) E S and A;_ (£, x) < j ^ A; + (i, x) 
— oo, j < k-(t,x) for all t such that (t,x) G 5 
+oo, j > h+(t,x) for all t such that (i, x) G 5. 

When |t(j, x) | ^ oo, (t(j, x), x) is an element of S. Given a fixed j, if t(j, xq) = +oo 
for some Xo > (resp. < 0) then t(j,x) = +oo for all x ^ Xo (resp. ^ Xo); 
same for — oo. For x = it is always the case that t(j, 0) G {— N, . . . , N}. 
Also, t(j, x) ^ t(j + 1, x) for j G {1, . . . , M - 1}, x G {—TV, . . . , jV}. For all 
j G {1, . . . , M - 1} there is x G {-N, ...,N} such that t(j, x) < £(j + 1, x). Also 
t(j, x + 1) = t(j, x) ± 1 for all j G {1, . . . , M}, x G {—TV, . . . ,N — 1} such that 
|t(j, x + l)|, |t(j, x)| 7^ oo. If x)| < oo and \t(j,x + l) \ = oo or \t(j,x — l)\ = oo 
then (t(j, x),x) G <96\ 
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To relate t(j,x) with k±(t,x) take t(0,x) = — oo, t(M + l,x) = +00 and 

fc_(t,x) = max{fc : < k < M,t(k,x) < t}, 

k+(t, x) = mm{k :0^k^M,t(k + l,x)> t}. 

Now for j G {1,...,M} take a = min{x : t(j,x) 7^ ±00}, b = maxji : 
t(j,x) 7^ ±00} and define yi = (t(j,x),x) for a ^ x ^ 6. Take £j = 
(2/a, e a+ i, y a+ i, . . . , e b , y b ), e j = (y^i, Then £ 1; . . . , £ M G C(5) and 

£1 -< ••• ^£ M - (A.5) 

Given the set l\ -< • • • -< £m, write -D(£j) = {xj ,xj + 1, . . . , —1, 0, 1, . . . , xj — 
1, xj} and y = (t,(x),x) for y G Then (tj(xf), xf) G <9±5 U d ± S. Take 

{i/(x), xj ^ X ^ x+ 
+00, x > xj and tj(xj) > or x < xj and tj(xj) > 
— 00, x > x+ and tj(xj) < or x < x~ and tj(xj) < 0. 

This completes the set of equivalences: 

fj <-)• J9 <->• (A, fc), <-> fcj- <-> t «-> 

For y = (t,x) G S 1 , take K(y) = {k~(y) + 1, . . . , k + (y)}. It follows from the 
definition of t(J, x) that 

for all (t,x) G S, j G {1,...,M}, j e K(t,x) t(j,x) = t. (A.6) 

For e G £(5), take fc-(e) = fc-(y) V k-(y') and fc+(e) = fc+(y) A k + (y'), where 
e = (y,y'). Also take K(e) = {k_(e) + 1, . . . , k+(e)} = K{y) n K(y')- Yet 
for £ = (yo,e 1 ,y 1 ,e 2 ,y 2 ,...,e n ,y n ) C S, take = K(y ) n • • • D K(y n ) = 

K( ei )n---nK(e n ). 

Claim A.3. Let £ C S and j G {1, . . . , M}. Then £ C ^- 1/ and on/y i/ j G K(f). 

Proof: Take £ C Write £ = (y a /, e a / + i, y a / +1 , . . . , y v ) and £j = 

(y a , e a+ i, y a +i, ■ ■ ■ ,yb) with a ^ a' < b' ^ b. By construction of £j, for all a ^ % ^ 6, 
t(j, jp.) = ti and, by ([Oil , j G # (j/i). Therefore j G if(a') fl • • • D = K{£). 

Now suppose j G K(£) for some £ C. S. Write £ = (y ,ei,y%, . . . ,y n ) and, for 
i = 0, ... ,n, yi — (ti,Xi). By definition j G K(yi) for all i, by ( 1A.6I) this implies 
t(j, Xj) = i j and by construction of £j we have y^ G ^ V i; it follows from this last 
fact that £ C £j, thus proving the claim. 
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Let J j = (qj-t, qj], j G {1, . . . , M}. Notice that 

jo n Ji' = when j ^ f. (A. 7) 

For y G S, take J{y) = \Jj E K(y) an d for £ C S" take J(£) = (JjeKY^) Given 
e = (y,y') G £(£), let J(e) = (g fe _( e ), < ?*,+(» ]■ Notice that J(e) = J(y) fl J(y'), 
since J(e) = U jeA - (e) J J '. It follows from flA~3ll - f[01) that 

\J(e)\ = V (e). 

By the above fact there is one translation T e that relates (0, 77(e)] with J(e). We 
associate the atoms and subintervals of these translated intervals according to the 
following rules. We write (ei,pi) ~ (e 2 ,p 2 ) when (e^T^pi) ~ (e 2) T~^p 2 ) and 
(ei, Ji) w (e 2 , J 2 ) when (e^T^Ji) ~ (e 2 ,r~ 1 J 2 ). 

Claim A. 4. Let e\ = (y ,yi),e 2 = {2/1, 2/2) with x < x\ < x 2 be two adjacent 
edges in £(S) and J\, J 2 intervals. Then (ei, J\) ~ (e 2 , J 2 ) i/ and onfo/ «/ J\ = J 2 
and Jx C J(y )nJ(y 1 )nJ(y 2 ). 

To prove the claim it suffices to show that (ei,pi) ~ (e 2 ,p 2 ) iff p\ = p 2 G J(ei) fl 
J(e 2 ), since T e : (0, 77(e)] — > J(e) is just a translation and J(ei) fl J(e 2 ) = </(t/o) H 

■/(yi) n J(t/ 2 ). 

So suppose (ei,pi) ~ (e 2 ,p 2 ). Write pi = T ei pi and p 2 = T e2 p 2 . Of course 
Pi G (0, 77(ei)], thus pi G J(ei); also p 2 G J(e 2 ). Writing T e more explicitly one 
gets T e p = p + qk_( e )- We have 4 cases to consider. Case 1: e\ = e^, e 2 = e^. 
In this case p\ = p 2 , k-(yo) ^ k-(yi) and k-(y 2 ) ^ k-(yi), thus fc_(ei) = A;_(e 2 ) 
and pi = pi + gfe_(ei) = p 2 + g/c_( e2 ) = P2, since by hypothesis pi = p 2 . Case 2: 
ei = e^, e 2 = e^. In this case p x = p2-r)( e O and M2/2) < MvOj thus fc 2 (e 2 ) = 
k-(yi). By (EED <?m?a) + 77(e^) = g^-i,^) + 77 (e£) = g fc ( tl)X1 _i) = g&_( ei ), so 
Pi = Pi + g fc _( ei ) = Pi + qk-fa) + Vicyj = P2 + 9k_(i/i) = P2- Cases 3 and 4 are 
similar. 

Conversely, suppose p\ = p 2 G J(ei) fl J(e 2 ). Consider Case 3, i.e., ei = e^, 
e 2 = e^; the other cases are similar. Takepi = T~ 1 p 1 = p\— qk_( ei ) = Pi—Qk-(yi) G 
(0, 77(d)] and p 2 = T~ l p 2 = p 2 - g Me2 ) = p 2 - g My2) = p 2 - qk^+i^+i) = 
P2-Px 1+ i, tl = P2-[Px u t 1 -i+v(e^ 1 )} =P2-gfc_( yi )-7/(e^) =Pi-7/(e^) G (0,7/(e 2 )]. 
Since p 2 > we have pi > 7/(e^). Therefore (ei,pi) ~ (e 2 ,p 2 ) and the claim holds. 

It will be convenient to work with a different representation of a broken line, 
which we dub translated broken line. The translated broken lines have a simpler 
representation that follows from Claim lA~4l Given a broken line 7 G -6(77) in 5, 
we define the object ip — ^(7) by 

^ = ((yo,ei,j/i, . . .,e n ,y n ), J) , 
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where (y , ex,yx,---, e n , y n ) = £(7) and J = T ei J x = T e , 2 J 2 = ■ ■ ■ = T Cn J n . Write 
~/(if>) for the unique broken line 7 such that if) = ipil)- For any if) of the above form 
we have ^(ip) = (yo, ex, J±, y±, . . . , e n , J n , y n ), where Jj = T~ X J. The translated 
broken lines have all the properties analogous to those already discussed for the 
broken lines. For if) = (£x, J), define £(if>) = £\ and w(if)) = \ J\ = w(j(if>)). 

Claim A. 5. Let £ = (y , ex, yi, . . . , e n , y n ) C S and J C (0, qju] ■ Then (ex, J) ~ 
(e„, J) if and only if J C J(£). 



The proof is short. Suppose (ei, J) ~ • - • ~ (e n ,J). Since (ex, J) ~ (e 2 , J), by 
Claim |AJJ we have J C J(y )n J{y x )n J(j/ 2 )- Also (e 2 , J) « (e 3 , J), thus J C J(y 3 ) 
as well. Similarly, for i = 2, . . . , n, (e^-i, J) ~ (e^, J) and by Claim IAT41 J C J(?/j). 
Therefore J C n" =0 J(^) = H™ =0 U ieK(2/i) P = U ie[n?=0 x( w )] J j = U ieXW J J = J(^); 
we have used ( 1A.7I) on the second equality. Conversely, if J C J(£) we have 
^ ^ J(Vi) fori = 0, . . . , n and by Claim IA.4I we have (ej_i, J) ~ (ej, J) for all 
z = 2, . . . , n, i.e., (ex, J) ~ • • • ~ (e n , J), which proves the claim. 

We define B(ff) = {ip(^f) : 7 G 5 (77)}. By definition of translated broken lines we 
have if) G B(ff) if and only if (ex, J) ~ ■ • • ~ (e n , J). It follows from Claim IA.5I 
that 

B(fj) = {iP = (I, J) :£C S, JC J(l)}. (A.8) 

As for the broken lines, given an £ C S, one can define the maximal translated 
broken line if>(£) in B(ff) that has trace £. It is clear that if)(£) = ip(j(£)) and 
7 (£) = j(if)(£)). Now by fTOl) one has if)(£) = (£, J(t)) and therefore 



w 



(£) = J2^jeK(,y (A.9) 



Let £ G C(£). By Claim [O] and Item [2] of Lemma Owe have j G if and 
only if £ = £j. Thus we have by ( 1A.5I) 

= <J y h 3 for any £ G C(5). (A.10) 
otherwise 



Now (J22U) follows from flATOl) and (jOl) . Finally, combining (jOl) and Claim [O] 

gives 

3 

and (I2.23P follows from (I2.2ip and the above equation. 

It remains to prove that, given any pair of sets {Ex -< £ 2 -< ■ ■ ■ -< £m} and 
{wx, • • • , % > 0}, equation (12.211) holds for some fj. For j = 0, 1, ... , M, let 
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qj = J2f=i w f an d let A = {0 = g < ?i < • ■ ■ < 9at}- For a; G {—AT, . . . , A/"}, 
j G {1, . . . , M} define t(J, x) from the {£j} as shown above. It then follows that 
t(j,x) has all the properties mentioned in the construction. From t(j, x), define 
k±{y),y G S and then K(y),y G S*. With A defined above and fc(-) we can 
recover (t, x) G 5* and from that we obtain the associated flow field fj. Now 
for such fj (12.2ip holds by the construction described in this appendix. 

A.5 Proof of ( E26p 

As in the proof of Lemma 12.11 take A(t,x) = {(t,x) : t ^ t + \x — x\}, A(V) = 
U y6 yA(y). Consider also A(t,x) = {(t,x) : t ^ t — \x — x\}. Then y G A(y') iff 
y' G It is clear that y m G so take n = mm{i : yi G There is 

y* G £ such that y n G -A (2/*) and thus y* G If n = 0, G S PI v4(y ) = 

{(xq, t ), (x + 1, to — 1), (x — 1, t — !)}• It is clear that any of these possibilities 
for imply t/ £ ^- So suppose n > 0. By construction x n _i — x n ± 1, so assume 
for simplicity x n „i = x n + 1. Now y n _i G" A(£), which means £ fl A(y n _ 1 ) = 0. 
So y* G A(y n )\A(y n -i) and thus x* ^ x n , t* — t n — (x n — a;*). Notice that 
£ must eventually reach {(t,x) : t — (t n + 1) = (x n + 1) — x} to cross S. Let 
y' = (t', x') be the point of the first time it happens, that is, the one with smallest 
x'. Since £ fl A(y n -i) = 0, we must have t' ^ t n + 1, so a;' ^ x n + 1. Thus 
t' — t* = (t 7 — t^) + (t^ — t^) ^ 1 + x n — x* ^ x' — x* and the equality holds if and 
only if t' — t n = 1, x' — x n = 1. But it must be the case that equality holds because 
( I2.19P implies \t' —t*\ ^ \x' — x*|. So y' — (t n + 1, x n + 1) G £. Now we only need 
to observe that £ must connect y* to j/' through y n . Indeed, if t(x n ) > t n we would 
have t(x n ) — t(x*) > x n — x* and if t(x n ) < t n we would have t(x') — t(x n ) > x' — x n 
and either of them is absurd because of (12.191) . 
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1.1 commutativity. The one-dimensional lattice with the particles' starting 
positions and the envelope sequences. Particles are represented by black 
balls below the axis and envelopes are represented by arrows (jump 
envelopes) or crosses (sleep envelopes) just above the axis. No matter 
in which order the evolution is performed, it will end up burning the 
envelopes below the bold line 

1.2 monotonicity. First, inserting a couple of sleep envelopes has decreased 
the final profile of burned envelopes. Second, starting with some par- 
ticles passive instead of active (squares instead of balls) has caused an 
analogous effect. Finally, even starting with quite a few extra sleep 
envelopes and some passive particles, the profile does not necessarily 
decrease 

2.1 a hexagonal domain 

2.2 construction of the brick diagram. In this example S is a 3 x 3 rectan- 
gular domain. There are 6 horizontal lines with a total of 24 intervals 
forming 9 bricks. The diagram is divided into 15 strips 

2.3 broken lines that cross S and calculation of the weight of a given broken 
trace. On the left we have a random configuration of broken lines on 
the 3x3 domain S and on the right we have the brick diagram from 
which they arose. The 3rd strip in the brick diagram is highlighted and 
the corresponding broken line appears in bold on the left. Also, for a 
broken trace £ that starts at the center of S, goes northeast, and then 
exits S going northeast again, we have determined which broken lines in 
C(S) contain £; they are the 11th and 12th ones and they also appear 
in bold. In this case w(£) equals the sum of the widths of the 11th and 
12th strips 
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construction of the maximal path. From £ one constructs the flow field 
and by following algorithm (|2.27p one gets an optimal path. The theory 
developed for broken lines and Theorem 12.11 guarantee this is indeed 
optimal. The algorithm forbids the path to cross any broken line; for 
that purpose it suffices to require that it does not 'cross' the flow field. |49 



Bibliography 



[AD95] D. Aldous and P. DlACONlS, Hammersley's interacting particle 
process and longest increasing subsequences, Probab. Theory Related 
Fields, 103 (1995), pp. 199-213. 

[And82] E. D. Andjel, Invariant measures for the zero range processes, 
Ann. Probab., 10 (1982), pp. 525-547. 

[AS89a] T. Arak and D. Surgailis, Markov fields with polygonal realiza- 
tions, Probab. Theory Related Fields, 80 (1989), pp. 543-579. 

[AS89b] , On polygonal Markov fields, in Stochastic methods in mathe- 
matics and physics, World Sci. Publ., 1989, pp. 302-309. (Karpacz, 
1988). 

[AS89c] , Polygonal fields: a new class of Markov fields on the plane, in 

Stochastic differential systems, vol. 126 of Lecture Notes in Control 
and Inform. Sci., Springer, Berlin, 1989, pp. 293-316. (Bad Honnef, 
1988). 

[AS89d] , Polygonal Markov random fields, Soviet Math. Dokl., 38 

(1989), pp. 284-287. 

[AS90] , Markov random graphs and polygonal fields with Y-shaped 

nodes, Probability theory and mathematical statistics, 1 (1990), 
pp. 57-67. 

[AS91] , Consistent polygonal fields, Probab. Theory Related Fields, 89 

(1991), pp. 319-346. 

[BTW87] P. Bak, C. Tang, and K. Wiesenfeld, Self- organized criticality: 
An explanation of the 1/ / noise, Phys. Rev. Lett., 59 (1987), p. 381. 

[BTW88] , Self- organized criticality, Phys. Rev. A, 38 (1988), pp. 364-374. 



73 



74 



BIBLIOGRAPHY 



[Dha99] D. Dhar, The abelian sandpile and related models, Physica A, 263 
(1999), p. 4. 

[DF91] P. Diaconis and W. Fulton, A growth model, a game, an alge- 
bra, Lagrange inversion, and characteristic classes, Rend. Sem. Mat. 
Univ. Politec. Torino, 49 (1991), pp. 95-119. (Commutative algebra 
and algebraic geometry II - Turin, 1990). 

[Dic02] R. Dickman, N onequilibrium phase transitions in epidemics and 
sandpiles, Physica A, 306 (2002), p. 90. 

[Dic08] , Restricted sandpile revisited, to appear in Phys. Rev. E, (2008). 

ArXiv:cond-mat/0601415| 

[DAM+01] R. Dickman, M. Alava, M. A. Munoz, J. Peltola, 
A. Vespignani, and S. Zapperi, Critical behavior of a one- 
dimensional fixed-energy stochastic sandpile, Phys. Rev. E, 64 
(2001), p. 56104. 

[DMVZ00] R. Dickman, M. A. Munoz, A. Vespignani, and S. Zapperi, 
Paths to self- organized criticality, Braz. J. Phys., 30 (2000), p. 27. 

[DTdO02] R. Dickman, T. Tome, and M. J. de Oliveira, Sandpiles with 
height restrictions, Phys. Rev. E, 66 (2002), p. 16111. 

[Eri96] K. Eriksson, Chip firing games on mutating graphs, SIAM J. Dis- 

crete Math., 9 (1996), pp. 118-128. 

[dFLdSHOO] J. E. de Freitas, L. S. Lucena, L. R. da Silva, and H. J. 

Hilhorst, Critical behavior of a two-species reaction- diffusion prob- 
lem, Phys. Rev. E, 61 (2000), pp. 6330-6336. 

[FMLOla] U. L. Fulco, D. N. Messias, and M. L. Lyra, Critical behavior 
of a one- dimensional diffusive epidemic process, Phys. Rev. E, 63 
(2001), p. 066118. 

[FMLOlb] , Critical dynamics of a Id diffusive epidemic propagation pro- 
cess, Physica A, 295 (2001), p. 49. 

[Gri95] G. Grinstein, Scale Invariance, Interfaces and N onequilibrium 
Dynamics, vol. 344 of NATO Advanced Study Institute, Series B: 
Physics, Plenum Press, New York, 1995. 



[Har78] T. E. Harris, Additive set-valued Markov processes and graphical 
methods, Ann. Probability, 6 (1978), pp. 355-378. 



BIBLIOGRAPHY 



75 



[HS] 

[HLM+08] 

[JanOl] 

[JPS98] 

[JohOO] 

[KS03] 

[KS05] 

[KS06] 

[KS08] 

[KSS89] 

[Lev04] 
[LP07] 

[Liib02a] 

[Liib02b] 

[LH03a] 



C. E. Hoffman and V. Sidoravicius. private communication. 

A. E. Holroyd, L. Levine, K. Meszaros, Y. Peres, 
J. PROPP, AND D. B. WILSON, Chip-firing and rotor-routing on 
directed graphs. arXiv:0801.3306V 2. 2008. 

H.-K. Janssen, Comment on "critical behavior of a two-species 
reaction- diffusion problem", Phys. Rev. E, 64 (2001), p. 058101. 

W. Jockusch, J. Propp, and P. Shor, Random domino tilings 



and the arctic circle theorem. |ArXiv:math/9801068, 1998. 



K. Johansson, Shape fluctuations and random matrices, Comm. 
Math. Phys., 209 (2000), pp. 437-476. 

H. Kesten and V. Sidoravicius, Branching random walk with 
catalysts, Electron. J. Probab., 8 (2003), p. 51. 

, The spread of a rumor or infection in a moving population, 



Ann. Probab., 33 (2005), pp. 2402-2462. 

, A phase transition in a model for the spread of an infection., 



Illinois J. Math., 50 (2006), pp. 547-634. 

, A shape theorem for the spread of an infection, to appear in 



Annals of Mathematics, 2008. 

R. Kree, B. Schaub, and B. Schmittmann, Effects of pollution 
on critical population dynamics, Phys. Rev. A, 39 (1989), pp. 2214- 
2221. 



L. LEVINE, The rotor-router model. arXiv:math/0409407vl, 2004. 



L. Levine and Y. Peres, Scaling limits for internal aggregation 
models with multiple sources. arXiv:0712.3378vl, 2007. 

S. Lubeck, Scaling behavior of the conserved transfer threshold pro- 
cess, Phys. Rev. E, 66 (2002), p. 046114. 

, Scaling behavior of the order parameter and its conjugated field 



in an absorbing phase transition around the upper critical dimension, 
Phys. Rev. E, 65 (2002), p. 046150. 

S. Lubeck and P. C. Heger, Universal finite-size scaling behavior 
and universal dynamical scaling behavior of absorbing phase transi- 
tions with a conserved field, Phys. Rev. E, 68 (2003), p. 056102. 



76 



BIBLIOGRAPHY 



[LH03b] , Universal scaling behavior at the upper critical dimension of 

n on equilibrium continuous phase transitions, Phys. Rev. Lett., 90 
(2003), p. 230601. 

[Man90] S. S. Manna, Large-scale simulation of avalanche cluster distribu- 
tion in sand pile model, J. Stat. Phys., 59 (1990), pp. 509-521. 

[Man91] , Two-state model of self-organized criticality, J. Phys. A: Math. 

Gen., 24 (1991), pp. L363-L369. 

[O'COO] N. O'Connell, Directed percolation and tandem 
queues, HP Labs Technical Reports HPL-BRLMS- 
2000-28, Hewlett-Packard, 2000. (downloaded from 



http://www.hpl.hp.com/techreports/2000/HPL-BRIMS-2000-28.html 



on February, 2008). 

[OvWLH89] K. Oerding, F. van Wijland, J. P. Leroy, and H. Hilhorst, 

Fluctuation-induced first-order transition in a nonequilibrium steady 
state, J. Stat. Phys., 99 (1989), pp. 1365-1395. 

[PW98] J. G. Propp and D. B. Wilson, How to get a perfectly random 
sample from a generic Markov chain and generate a random spanning 
tree of a directed graph, J. Algorithms, 27 (1998), pp. 170-217. (7th 
Annual ACM-SIAM Symposium on Discrete Algorithms - Atlanta, 
GA, 1996). 

[RPSV00] M. Rossi, R. Pastor- Satorras, and A. Vespignani, Uni- 
versality class of absorbing phase transitions with a conserved field, 
Phys. Rev. Lett., 85 (2000), pp. 1803-1806. 

[Ros] H. Rost. unpublished. 

[Ros81] , Nonequilibrium behaviour of a many particle process: density 

profile and local equilibria, Z. Wahrsch. Verw. Gebiete, 58 (1981), 
pp. 41-53. 

[Sep98] T. Seppalainen, Coupling the totally asymmetric simple exclusion 
process with a moving interface, Markov Process and Related Fields, 
4 (1998), pp. 593-628. I Escola Brasileira de Probabilidade (Rio de 
Janeiro, 1997). 



[She] 



E. Shellef. private communication. 



BIBLIOGRAPHY 



77 



[SSV] V. Sidoravicius, D. Surgailis, and M. E. Vares, Discrete 

geometric broken line process. Preprint. 

[SSV99] , Poisson broken lines process and its application to Bernoulli 

first passage percolation., Acta Appl. Math., 58 (1999), pp. 311-325. 
(Limit theorems of probability theory - Vilnius, 1999). 

[vWOH98] F. van Wijland, K. Oerding, and H. J. Hilhorst, Wil- 
son renormalization of a reaction- diffusion process, Physica A, 251 
(1998), pp. 179-201. 



